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1 Introduction 



The main goals of the Bhabha working group are to make an inventory of all the available 
Monte Carlo (MC) event generators for small-angle (SABH) and large-angle (LABH) Bhabha 
processes at LEP1 and LEP2, and to improve our understanding of their theoretical uncertain- 
ties through systematic comparisons of the MC event generators (developed independently) 
among themselves and with other non-MC programs. The presented activity is of course an 
obvious continuation of the previous workshops on LEP1 physics [1,2]. In the beginning of 
the present workshop the theoretical uncertainty at LEP1 for the SABH process was typically 
estimated as 0.16%, and for the LABH process was estimated at 0.2% level at Z peak and 1% 
on the wings of the Z resonance. There were no estimates specific to LEP2. 

We shall concentrate on the comparison of all the presently available theoretical calculations 
(published and unpublished). This will be done for several kinds of event selection (ES), defined 
as a set of experimental cuts and apparatus acceptances, starting from ES's unrealistic, but 
useful for some studies oriented towards the QED matrix element, and ending on ES's very 
close to the experimental ones. 

Let us add a few comments to clarify our priorities and to set the proper perspective for our 
work. In spite of the considerable effort of several theoretical groups, at present the theoretical 
error on the small-angle Bhabha cross section dominates the luminosity error at LEP1. This 
inhibits from taking full advantage of the high experimental precision of the final LEP1 data 
for precision tests of the Standard Model. As a consequence, the reduction of the theoretical 
error in the SABH process at LEP1 is the biggest challenge, and was the main objective of 
our working group. The precision requirements of LEP2 are lower than those of LEP1. The 
total cross section of W pair production will be measured with 1.0% to 0.5% precision at best, 
so it is sufficient to keep the theoretical uncertainty of the SABH process at the 0.25% level. 
Furthermore, at LEP2 the detectors and experimental techniques for measuring the SABH 
process are almost the same as for LEP1 1 . Radiative corrections to the SABH cross section 
depend on the center of mass energy, but smoothly; moreover, in the small-angle regime the 
center of mass energy is not so important from the point of view of the physics involved: we are 
always faced with a ^-channel photon-exchange dominated process; hence, improving the small- 
angle Bhabha generators for LEP1 is generally a sufficient condition for improving them also 
for LEP2. The only subtle point concerns the error estimate: a 0.1% error at LEP1 does not 
guarantee such a small error also at LEP2, so that an additional analysis has to be performed. 
For the LABH process, the final LEP1 data analysis requires a theoretical uncertainty of the 
codes used to be at the 0.5% level. The LABH process at LEP2 is not of major interest, 
and we think that a precision of the order of 2% is enough. Nevertheless, the physics of the 
LABH process at LEP2 is significantly different from LEP1 (different Feynman diagrams rise to 
importance), so performing additional study for the LABH process at LEP2 is a new nontrivial 

1 Actually, the main difference is that, due to machine background radiation, the internal part of luminosity 
detectors may be obscured by special masks. We shall discuss the impact of such modification on the theoretical 
errors. This aspect was brought to our attention by B. Bloch-Devaux during our WG meeting in January 95. 
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work 2 . 



In view of the above, our strategy was to do all the work for the SABH and the LABH 
processes first for LEP1 experimental conditions, and to supplement it with all necessary 
work/discussion which would assure control of the precision at the level sufficient for LEP2 
experiments. This practically means that all the numerical comparisons were done for LEP1 
and repeated for LEP2, or, in rare cases, a convincing argument was given that it is not nec- 
essary (sometimes numerical results for LEP2 were obtained, but are not shown in full form 
because they were trivially identical to those for LEP1). 

We include in our report two main parts: one part on the SABH process and a second one on 
the LABH process, with the cases of LEP1 and LEP2 discussed in parallel. These two processes 
are governed by different physics (i.e. dominated by different Feynman diagrams). Also, the 
theoretical precision requirements in calculating SABH and LABH cross sections are different 
by a factor of five-ten. These two parts are followed by a section including short descriptions 
of all the involved Monte Carlo (MC) event generators or other codes, and a final section on 
conclusions and outlook. 



2 Small-angle Bhabha scattering 

Small-angle Bhabha (SABH) scattering is used at LEP1 and LEP2 to measure the accelerator 
luminosity. The LEP1 experiments have reached in 1993-1994 a systematic uncertainty of 
better than 0.10% in selecting luminosity Bhabha events, see Ref. [3] and Refs. [4,5]. 

On the theory side, QED calculations have still an uncertainty larger than 0.16% [6] in 
determining the Bhabha cross section in the detector acceptance, which is caused mostly by 
the non-existence of a Monte Carlo program including complete 0(a 2 ) next-to-leading terms. 
Actually, there exist 0(a 2 ) calculations with complete next-to-leading contributions [7,8] which 
claim a precision of the order of 0.1%, but they can not be used in a straightforward way, 
because they are not implemented in the Monte Carlo event generators. The size of the 0(a 2 ) 
contributions depends not only on the angular range covered by the detector and on the electron 
energy cut-off, but also on crucial experimental aspects, such as the sensitivity to soft photons 
or such as the electron cluster size. This means that the main interest is in the theoretical 
predictions for the Bhabha process, including as many higher order radiative corrections at it 
is necessary to reach a precision of 0.05%, in a form of a Monte Carlo event generator. 

Monte Carlo event generators are very powerful tools because they are able to provide a 
theoretical prediction, cross sections and any kind of distribution, for arbitrary ES's. However, 
event generators are difficult to construct and, what is even more serious, they are very difficult 
to test - one has to have at least two of them to compare with one another for a wide range of 

2 The radiative LABH process is an important background to other processes, like r-pair production, 
W + W~ — > eei/v, "new physics" like SUSY processes and so on, but a detailed analysis of these items goes 
beyond the aims of the present study. 
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ES's. 



For the SABH process, the task of comparing various Monte Carlo event generators was the 
main goaloi the Bhabha Working Group. There were only a few comparisons of independently 
developed Monte Carlo event generators for the SABH process in the past. A few examples can 
be found in Ref. [2]. However, we shall include in the comparisons results from non-Monte Carlo 
calculations, as well. They are usually limited to certain special (primitive) ES's. Nevertheless 
they provide additional valuable cross-checks. 

What shall we learn from these comparisons? The calculations from various Monte Carlo 
event generators will of course differ. The differences have to be understood. In a certain class 
of the comparisons, the underlying QED matrix element will be the same and in that case 
the differences will be only due to numerical effects. The results from two or more computer 
programs will differ due to rounding errors, programming bugs, numerical approximations. The 
difference measures uncertainties of this kind, and we say that we are determining the technical 
precision of the tested programs. One has to remember that the technical precision is dependent 
on the ES, and it is therefore absolutely necessary to use several at least semi-realistic, quite 
different, examples of ES's. In other cases, we shall compare Monte Carlo event generators 
which are based on different QED matrix elements. In this case, the difference between results 
will tell us typically about higher order effects which are not included in some of these event 
generators, or which are approximated differently in these programs. In this situation we shall 
talk about exploring the physical precision of the tested Monte Carlo event generators. Needless 
to say, the physical precision is the main goal, but one has to remember that without a technical 
precision of at least a factor of two better than the physical precision it is pointless to discuss 
the physical precision at all! 

Before we come to the actual comparisons of the programs, let us characterize various 
contributions/corrections to the SABH process. We shall also characterize briefly the various 
Monte Carlo event generators and non-Monte-Carlo calculations involved in the comparisons. 

If we remember that the SABH process was chosen for the luminosity measurement because 
it is calculable from first principles within quantum electrodynamics (QED), then it is natural 
to group corrections to the SABH process into pure-QED and non-QED corrections. The 
latter ones are due to s-channel Z-exchange, and the corrections induced by low energy strong 
interactions (QCD) through vacuum polarization and light quark pair production. Among the 
pure QED corrections, we may distinguish photonic (bremsstrahlung) corrections, related to 
multiple photon emission, and non-photonic corrections - for instance lepton pairs, leptonic 
vacuum polarization, multiperipheral diagrams. Numerically, the biggest ones are the photonic 
corrections and the vacuum polarization correction. They also contribute the most to the 
physical precision. Photonic corrections dominate completely the technical precision, due to 
the MC integration over the complicated multi-body phase space. The QED non-photonic 
corrections are small, but are difficult to calculate and quite uncertain (technical precision). 

For all the comparisons of the event generators it is crucial (especially for SABH) to under- 
stand the experimental ES. In the main comparisons we shall compare all the available event 
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generators for four types of ES's. However, the problem of the variation of the parameters in 
the ES is so important that we include also a separate subsection on this subject, in which, for 
a limited number of three event generators, we perform a detailed study of the dependence of 
the higher order corrections on all possible cut-offs involved in the real experiment. This will 
allow us to see all our work in the proper perspective from the point of view of the experimental 
analysis, and will also give us clear hints on the dependence of the higher order corrections on 
the fine details of the ES. This study will be limited to the SABH process. 

2.1 Sensitivity of LEP1 observables to luminosity 

The importance of the improvement of the theoretical luminosity error on the LEP1 results 
is shown in Table 1. The results of the lineshape parameter fits made with the theoretical 
luminosity error of 0.16% and 0.11% are given [9], corresponding to the reduction of error 
achieved during this workshop. A projection concerning a further reduction of the theoretical 
luminosity error to 0.06% is also given. The results of the four LEP1 experiments used as input 
to the fits, as well as the fitting procedure, are described in Ref. [10]. From the five parameter 
fit, only <r° is sensitive to the luminosity error. The decreased error in this variable causes 
a reduction of the errors of the derived parameters shown in the lower part of Table 1. As 
we see, the above improvement in the theoretical luminosity error influences significantly not 
only quantities like the "number of light neutrino's" N v , but also other LEP1 observables used 
routinely in the tests of the Standard Model. 





theoretical luminosity error 




0.16% 


0.11% 


0.06% 


mz [GeV] 


91.1884 ± 0.0022 


91.1884 ± 0.0022 


91.1884 ± 0.0022 


T z [GeV] 


2.4962 ± 0.0032 


2.4962 ± 0.0032 


2.4961 ± 0.0032 


< [nb] 


41.487 ±0.075 


41.487 ±0.057 


41.487 ± 0.044 


Ri 


20.788 ± 0.032 


20.787 ±0.032 


20.786 ± 0.032 


A ' 1 

FB 


0.0173 ±0.0012 


0.0173 ±0.0012 


0.0173 ±0012 


T had [GeV] 


1.7447 ± 0.0030 


1.7447 ± 0.0028 


1.7446 ± 0.0027 


r„ [MeV] 


83.93 ±0.13 


83.93 ±0.13 


83.93 ±0.12 


*S [nb] 


1.9957 ±0.0044 


1.9958 ± 0.0038 


1.9959 ± 0.0034 


Thad/Tz [%] 


69.90 ±0.089 


69.90 ±0.079 


69.89 ±0.072 


T u /T z [%] 


3.362 ± 0.0037 


3.362 ± 0.0032 


3.362 ± 0.0028 


T inv [MeV] 


499.9 ± 2.4 


499.9 ±2.1 


499.9 ± 1.9 


Finv/Fu [%] 


5.956 ±0.030 


5.956 ± 0.024 


5.956 ±0.020 


N v 


2.990 ±0.015 


2.990 ±0.013 


2.990 ±0.011 



Table 1: Line shape and asymmetry parameters from 5-parameter fits to the data of the four LEP1 
experiments, made with a theoretical luminosity error of 0.16%, 0.11% and 0.06% [9]. In the lower part of 
the Table also derived parameters are listed. 

At LEP2, the normalization of the total cross section for the WW production process enters 
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in a nontrivial way into tests of the W boson coupling constants. The precision requirements 
for the total cross section is limited by statistics of the WW process, and a luminosity error 
at the 0.25% level is sufficient (see the chapter "WW cross-sections and distributions", these 
proceedings). 

2.2 Higher order photonic corrections at LEP1 and LEP2 



Canonical coefficients 




Omin = 30 mrad 


Omin = 60 mrad 




LEP1 


LEP2 


LEP1 


LEP2 


O(aL) 




137xl0" 3 


152xl0" 3 


150xl0" 3 


165xl0" 3 


0(a) 


9 1 a 

2 7T 


2.3xl0" 3 


2.3xl0" 3 


2.3xl0" 3 


2.3xl0" 3 


0(a 2 L 2 ) 


Hf«) 2 


9.4xl0" 3 


llxlO" 3 


llxlO" 3 


14xl0" 3 


0(a 2 L) 


SL (SL^L) 

7T \7T / 


0.31xl0" 3 


0.35xl0" 3 


0.35xl0" 3 


0.38xl0" 3 


0(a 3 L 3 ) 


i (;«) s 


0.42 xlO" 3 


0.58xl0" 3 


0.57xl0" 3 


0.74xl0" 3 



Table 2: The canonical coefficients indicating the generic magnitude of various leading and subleading 
contributions up to third-order. The big-log L = ln(|f|/m 2 ) — 1 is calculated for 9 m in = 30 mrad and 
9 m i n = 60 mrad and for two values of the center of mass energy: at LEP1 (^/s = Mz), where the 
corresponding \t\ = (s/4)^ ifl are 1.86 and 7.53 GeV 2 , and at LEP2 energy (y/s = 200 GeV), where the 
corresponding \t\ are 9 and 36 GeV 2 , respectively. 

For the SABH process, the smallness of the electron mass "ruins" the normal perturbative 
expansion order in the following sense: for instance, the 0(a 2 ) QED contributions can be 
expanded into 0(a 2 L 2 ), 0(a 2 L) and pure non-log 0(a 2 ). The non-log 0(a 2 ) corrections 
are completely uninteresting, while the 0(a 3 L 3 ) corrections are as important as the 0(a 2 L) 
corrections. Here L = ln(|£|/ra 2 ) is the so-called big-log in the leading-logarithmic (LL) ap- 
proximation, where t is the momentum transfer in the ^-channel (of the order of 1 GeV). This 
phenomenon is illustrated in Tab. 2. From this table, it is clear that for a precision of the order 
of 0.25% (for calorimetric ES's) it is enough to include the 0(a 1 L), 0(a 1 ) and 0(a 2 L 2 ). For 
a precision of the order of 0.1% or better, one has to add 0(a 3 L 3 ) and 0(a 2 L). These "scale 
coefficients" have to be kept in mind when discussing various QED calculations/programs. As 
we shall see, the higher order effects seen in the numerical results presented in the next sections 
generally conform to the above scale coefficients. 

Table 2 demonstrates also the "scaling laws" for various QED corrections between LEP1 
(Z peak) and LEP2 energies. If the angular range is kept the same, then ^-channel transfer is 
proportional to s = 4:E 2 eam . Actually, at LEP2 experiments the luminosity measurement will 
rely more on the SABH process at larger angles, above 3°, and this is why we also included in 
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the table another two columns for this angular range. As we see, photonic corrections do not 
change very much due to the increase of y^s from Z-peak energy to LEP2 energy (200 GeV) and 
due to going to twice larger angles. Actually, the change in canonical coefficients is negligible. 
One has only to pay attention to the 0(a 3 L 3 ) corrections, which in the worst case increase by 
a factor 1.75 (however, as we shall see they are under good control). 

One has to remember that, as it was shown explicitly in ref. [11], the radiative corrections to 
the SABH process with the typical "double tag" detection are proportional to ln((# max /# m i n ) — 
1), i.e. they are bigger for "narrower" angular acceptance and smaller for "wider" angular ac- 
ceptance. This has to be remembered, because at LEP2 in some experiments the angular range 
might be "narrowed" by placing masks in front of the SABH detectors in order to eliminate 
machine background radiation. We conclude that the change for "narrower" angular acceptance 
is more dangerous from the point of view of the increase of the pure photonic corrections, and 
we shall address this problem with a separate numerical exercise. 

In ref. [11] it was also shown (numerically), using an O(a) calculation, that for the purpose 
of the SABH process below 6° we may neglect the real and virtual QED interference contri- 
butions between photon emission from the electron and positron lines, the so called "up-down 
interference". In the numerical example in ref. [11] it was shown that, for the angular range 
3.0° — 4.24°, the "up-down interference" is below 0.015%. It is even smaller for smaller an- 
gles. It means that it is negligible for all practical purposes in the luminosity measurements. 
This phenomenon was also discussed in ref. [12] beyond O(a) in the framework of the eikonal 
approximation. 

2.3 Light pairs and other small contributions 

To calculate pair corrections to the SABH two approaches have been used. (1) The first one 
is based on direct analysis of Feynman graphs and analytical extraction of graphs and terms 
contributing to the SABH within the (9(0.1%) accuracy. Both leading and next-to-leading 
terms are considered. (2) The other method uses the LL approximation to find the dominant 
pair contributions to SABH and to discard the negligible ones. Having isolated the dominant 
mechanism, an actual MC program for this particular mechanism is constructed. 

(1) The dominant pair production corrections (enhanced by factors of L 2 and L) arise from 
kinematical configurations where one (or both) of the produced leptons is almost collinear with 
the incoming or outgoing e ± . These contributions have been calculated analytically [13,14]. 
The analytical calculation [7, 8, 13-15] of the real hard pair production cross-section within 
logarithmic accuracy takes into account the contributions of the collinear and semi- collinear 
kinematical regions. All possible mechanisms for pair creation (Singlet and Non- Singlet), as 
well as the identity of the particles in the final state, are taken into account 3 . In the case of 

3 Here we have taken into account only e + e~ pair production. An estimate of the muon pair contribution gives 
less than 0.05% since \a(Q 2 /m 2 ) ~ 3MQ 2 / m D- Contributions of pion and tau-lepton pairs give still smaller 
corrections. Therefore, within the 0.1% accuracy, one may omit any pair production contribution except the 
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Figure 1: The Feynman diagrams giving logarithmically enhanced contributions in the kinematical region 
where the created pair goes along the electron direction. The signs represent the Fermi-Dirac statistics of 
the interchanged fermions. 



Channel 


e e 




TT 


cc 


uu, dd, ss 


total 


a (nb) 


0.006 


0.006 


0.0008 


0.0005 


0.0011 


0.0144 



Table 3: Double Tag cross sections for fermion pair production from multiperipheral graphs. y/s = 
91.2 GeV, 30 mrad < 6 e +,6 e - < 60 mrad. For u,d,s quarks W 77 > 4 GeV. The uncorrected Born cross 
section a\B orn is 104 nb. 



SABH only a part of the total 36 Feynman diagrams are relevant, i.e. the scattering diagrams 4 
shown in Fig. 1. 

The analytical formulae for virtual, soft, hard and total pair production contributions can be 
found in [7,8,15]. Numerical results for the pair contribution cross sections based on these 
formulae can be obtained by using the code MLLBHA (see below for a description of the code). 
The leading term can be described by the electron structure function D e e [[x) [16-22]. Numerical 
results can be found in Refs. [2,7,8,15,23]. The contribution to SABH of the process of pair 
production accompanied by photon emission when both, pair and photons, may be real and 
virtual has also been analyzed and the relevant analytical formulae are given in [2,23]. 

With the help of a Monte Carlo generator [24,25], a dedicated study has been done for the 
contribution of the multiperipheral graphs, Fig. 1 (5-8), being for many kinematical setups the 
dominant mechanism of pair production. The total cross sections for the production of fermion 
pairs as detailed in Table 3 were obtained. The total contribution from the multiperipheral 
graphs is then estimated to be 1.4 X 10~ 4 <TB orn , with a relative error (from MC statistics) of 

e + e~ one. 

4 It can be verified [7, 8, 15] that the interference between the amplitudes describing the production of 
pairs moving in the electron direction and the positron one cancels. This is known as up-down (interference) 
cancelation. 
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0.3 


0.4 


0.5 


0.6 


0.7 


.3/175GeV 


.7/175GeV 


10 4 <T NN I CTBorn 
^ A(T Wwl a Born 
lO^aj^y/ aBorn 


-3.619 
-2.748 
-2.142 


-3.655 
-2.798 
-2.191 


-3.707 
-2.883 
-2.264 


-3.807 
-3.175 
-2.478 


-4.191 
-3.771 
-3.064 


-4.185 
-3.177 
-2.489 


-4.884 
-4.405 
-3.603 



Table 4: LL Non-Singlet ee pair correction to SABH, SiCAL angular cuts WW: 1.5° - 3.15°, NN: 1.61° - 
2.8°, in 10 4 X Born units, ^/s = 91.1888 GeV (175 GeV for last two entries), z m { n = s' m i n /s. 



±50%. This correction, which still does not take into account a further reduction factor of 
~ 20 coming from a cut on the acoplanarity angle of the detected e ± , is thus negligible for 
SABH [26]. 

(2) The LL calculation of photonic corrections to SABH of Ref. [27] has been extended 
to pair corrections in [28]. Analytical formulae for arbitrary asymmetric angular cuts, for 
both Singlet and Non-Singlet corrections have been given in [28] 5 . These formulae, based 
on [22], include both pairs and photons up to the exponentiated second or third order. The 
semianalytical program BHPAIR based on this calculation has been written [28]. Numerical 
results for the LCAL type angular cuts have been given in [28]. For the SiCAL type angular 
cuts the Singlet contribution is negligible (below 5 X 10~ 5 <TB orn ) and the Non-Singlet contribution 
(up to third order with exponentiation) is calculated in Table 4, also for the LEP2 energies. 
The strong dependence of the result on angular cuts (WW, NN or NW) may indicate significant 
effects due to more realistic ES's. This can only be analyzed with the MC simulation. Such 
a MC program has already been constructed [29]. This program, being an extension of the 
BHLUMI MC code [30], is based on the extension of the YFS resummation of soft photons to 
the resummation of infrared and collinear pairs, cf. [31]. Preliminary results [29] show that a 
calorimetric ES reduces further the pair correction of Table 4. 

To summarize, numerical values of pair corrections as given in [7,8,15], [2,23] and Table 4 
agree within 4 X 10~ 4 <TB orn for the NN and WW cuts. The total contribution from pairs and 
multiperipheral diagrams for the energy cut in the experimentaly interesting range 0.3 < x c < 
0.7 is also at most 4 X 10~ 4 <TB orn . With the help of a MC simulation of a realistic ES, one 
should be able to control the pair contribution with an accuracy of 3 X 10~ 4 <TB orn , or better. A 
similar conclusion is to be expected also for the LEP2 energies. 

2.4 Vacuum polarization 

Vacuum polarization contributes about 5.3% and 4%, respectively, to the e + e~ cross-section 
in the angular region of the first and second generation of the luminosity detectors at LEP 
[3,32]. The leptonic part of this contribution is known with excellent precision. The quark 

5 Extending further the analysis of Ref. [28], with the help of the 'parton-like' picture together with appro- 
priate choices of structure functions and hard scattering cross-sections, one can calculate the other pair creation 
mechanisms, including the multiperipheral one, as well as other leptonic backgrounds to SABH resulting from 
the 'charge blindness' of the detectors. This analysis will appear elsewhere [29]. 
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part, however, is more difficult since the quark masses are not unambiguously defined and 
perturbative QCD cannot be used for reliable calculations [33-35]. Therefore this part is 
calculated using a dispersion integral of Rhad 

a(e + e~ -^-hadrons) 

Rhad = , , _ , _x (1) 

a{e + e — >■ fi+fi ) 

measured experimentally. 






\t\ 




(a) 




(b) 




(rad) 


(GeV 2 ) 


Ref. [36] 


Ref. [33] 




Ref. [34] 




.020 


.83 


-.00345 


-.00340 ± .00008(2.£ 


%) 


-.00339 ± .00013(3.9%) 


.00002 


.030 


1.87 


-.00505 


-.00494 ± .00014(2i 


%) 


-.00493 ± .00020(4.1%) 


.00001 


.040 


3.33 


-.00629 


-.00612 ± .00019(3.1%) 


-.00613 ± .00025(4.1%) 


-.00001 


.050 


5.20 


-.00729 


-.00711 ± .00024(3.4%) 


-.00714 ± .00030(4.2%) 


-.00003 


.060 


7.48 


-.00812 


-.00795 ± .00027(3.£ 


%) 


-.00801 ± .00034(4.3%) 


-.00006 


.070 


10.18 


-.00889 


-.00869 ± .00030(3.£ 


%) 


-.00876 ± .00038(4.4%) 


-.00006 


.080 


13.30 


-.00963 


-.00936 ± .00033(3.£ 


%) 


-.00941 ± .00040(4.3%) 


-.00005 


.090 


16.83 


-.01029 


-.00997 ± .00035(3.£ 


%) 


-.01000 ± .00043(4.3%) 


-.00003 


.100 


20.77 


-.01089 


-.01052 ± .00037(3.£ 


%) 


-.01058 ± .00045(4.3%) 


-.00006 


.110 


25.13 


-.01144 


-.01103 ± .00039(3.£ 


%) 


-.01110 ± .00047(4.2%) 


-.00007 


.120 


29.90 


-.01194 


-.01150 ± .00040(3.£ 


%) 


-.01157 ± .00049(4.2%) 


-.00008 


.130 


35.08 


-.01241 


-.01193 ± .00042(3.£ 


%) 


-.01201 ± .00050(4.2%) 


-.00008 



Table 5: The hadronic part of the vacuum polarization contribution to the small-angle Bhabha scattering 
as a function of the scattering angle (and corresponding momentum transfer t). In column 4 and 5 also 
the ratio of the error to the value of the hadronic contribution is given in brackets. The last column gives 
the difference between the results of Refs. [34] and [33]. 



Recently, several reevaluations of the hadronic contribution to the QED vacuum polarization 
have been performed, mainly to determine the effective QED coupling a(m 2 z ) [33,37-42] and 
the anomalous magnetic moment (g-2) of the leptons [33]. At the same time, the vacuum 
polarization contribution to the small-angle Bhabha scattering has been recalculated [33,34]. 
Table 5 compares the results of these two calculations of the hadronic contribution in the 
angular region of small-angle Bhabha scattering used at LEP for the luminosity measurements. 
They are in excellent agreement, as is evident from the very small differences listed in the last 
column. In brackets, the error is given as a percentage of the total hadronic contribution. We 
see that the error of Ref. [33] varies between 63% and 83% of that of Ref. [34] in the angular 
region presented here. Numbers have been obtained with the help of FORTRAN routines 
HADR5 [33] and REPI [34] available from the authors. Finally the values of the previously 
used hadronic contribution from Ref. [36] are also shown. 

Fig. 2 from Ref. [34] shows the contribution of different energy regions of R to the value of 
the hadronic contribution and its error while the Fig. 3 from Ref. [33] shows the uncertainty 
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Figure 2: Relative contributions to Aa(f = —1.424 GeV 2 ) in magnitude and uncertainty from the Ref. [34]. 
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Figure 3: Relative uncertainty in percent of the hadronic vacuum polarization contribution as a function 
of the momentum transfer in the small-angle Bhabha scattering calculation from the Ref. [33]. 
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(rad) 


Ref. [36] 


Ref. [33] 


Ref. [34] 


hadronic / 07 \ 
total l /0 i 


.020 


-.01590 


-.01585 ± .00008 


-.01583 ± .00013 


21 


.030 


-.01877 


-.01866 ± .00014 


-.01865 ± .00020 


26 


.040 


-.02095 


-.02078 ± .00019 


-.02079 ± .00025 


30 


.050 


-.02271 


-.02252 ± .00024 


-.02255 ± .00030 


32 


.060 


-.02418 


-.02400 ± .00027 


-.02406 ± .00034 


33 


.070 


-.02551 


-.02531 ± .00030 


-.02537 ± .00038 


35 


.080 


-.02674 


-.02647 ± .00033 


-.02652 ± .00040 


36 


.090 


-.02785 


-.02753 ± .00035 


-.02756 ± .00043 


36 


.100 


-.02886 


-.02849 ± .00037 


-.02855 ± .00045 


37 


.110 


-.02979 


-.02938 ± .00039 


-.02945 ± .00047 


38 


.120 


-.03064 


-.03020 ± .00040 


-.03028 ± .00049 


38 


.130 


-.03144 


-.03096 ± .00042 


-.03104 ± .00050 


39 



Table 6: The vacuum polarization contribution to the small-angle Bhabha scattering as a function of the 
scattering angle. The last column gives the ratio of the hadronic part to the total vacuum polarization 
contribution. 



Generation 


typical (rad) 


Ref. [33] 


Ref. [34] 


first 
second 


.060 
.030 


.0003 
.0005 


.0004 
.0007 



Table 7: Summary of the uncertainty of the vacuum polarization calculation for the first and second 
generation of the luminosity detectors of LEP according to Ref. [33,34]. 

of the hadronic vacuum polarization contribution to the calculation of the small-angle Bhabha 
scattering as a function of the momentum transfer. 

The total vacuum polarization contribution is obtained as sum of the leptonic contribution 
and the hadronic one. It is shown in Table 6. The contribution of the vacuum polarization error 
to the total error of the luminosity measurement is about twice the error given in the Table 6. 
The typical angular region of the first and second generation of the LEP luminosity detectors 
is 60 and 30 mrad, respectively [3]. The contribution of the vacuum polarization error to the 
luminosity calculation for the LEP detector is given in Table 7. 

The vacuum polarization correction and its uncertainty are smaller for the lower angles 
covered by the second generation of luminosity detectors. 

In conclusion, the error of the hadronic contribution of Ref. [34] makes a negligible contri- 
bution to the total error of the calculation of the small-angle Bhabha scattering. The error of 
Ref. [33] is even smaller. Thus the error of Ref. [34] can be considered as a conservative one. 
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2.5 Brief characteristics of the programs/calculations 



Here we will very briefly summarize the basic features of the codes involved in the SABH 
comparisons. The only aim of the following is to just settle the frame, and not to give an 
exhaustive description of the codes, which can be found in the original literature and/or in the 
dedicated write-up's at the end of the present report. 

BHAGEN95 [43] - It is a Monte Carlo integrator for both small- and large-angle Bhabha 
scattering. It is a structure function based program for all orders resummation, including 
complete photonic O(a) and leading logarithmic 0(a 2 L 2 ) corrections in all channels. 

BHLUMI [44] - Full scale Monte Carlo event generator for small-angle Bhabha scattering. It 
includes multiphoton radiation in the framework of YFS exclusive exponentiation. Its matrix 
element includes complete O(a) and 0(a 2 L 2 ). The program provides the full event in terms 
of particle flavors and their four-momenta with an arbitrary number of radiative photons. 

LUMLOG - It is a Monte Carlo event generator for SABH (part of BHLUMI, see [44]). Photonic 
corrections are treated at the leading logarithmic level at the strictly collinear and inclusive way. 
Structure functions exponentiated up to G(a 3 L 3 ) are included (and without exponentiation up 
to 0(a 2 L 2 )). 

NLLBHA [2,23] - It is the FORTRAN translation of a fully analytical up to 0(a 2 ) calculation, 
including all the next-to-leading corrections. It is also able to provide 0(a 3 L 3 ) photonic cor- 
rections and light pair corrections including simultaneous photon and light pair emission. Not 
an event generator. 

OLDBIS - Classical Monte Carlo event generator for SABH from PETRA times [45] (the 
modernized version is incorporated in the BHLUMI set [44]). It includes photonic corrections 
at the exact 0(a). 

OLDBIS+LUMLOG - It is the well known "tandem" developed in order to take into account 
higher order corrections (LUMLOG) on top of the exact 0(a) result (OLDBIS). The matching 
between 0(a) and higher orders is realized in an additive form. 

SABSPV [46] - It is a new Monte Carlo integrator, designed for small-angle Bhabha scattering. 
It is based on a proper matching of the exact 0(a) cross section for ^-channel photon exchange 
and of the leading logarithmic results in the structure function approach. The matching is 
performed in a factorized form, in order to preserve the classical limit. 

2.6 Experimental event selection and theory uncertainty in lumi- 
nosity measurements 

In this section we discuss the interplay between experimental selection and higher-order radia- 
tive corrections. All numerical examples are for LEP1 at Z peak energy. The discussion of 
the results is generally limited to LEP1 but using "scaling rules" from the introduction one 
may easily extend it LEP2. In particular one has to remember that third order LL corrections 
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have the strongest energy dependence, and going from the Z-peak to the highest LEP2 energy 
introduces in them a factor of almost two. 

In this subsection three different event generators are used: i) a generator based on a com- 
plete first-order calculation OLDBIS [6,30], which has at most one photon radiated; it includes 
O(a) and O(aL); ii) a generator based on a leading-logarithmic third-order exponentiated 
calculation LUMLOG [6,30]; it includes O(aL), 0(a 2 L 2 ), 0(a 2 L 2 ) in strictly collinear ap- 
proximation; the 4-momenta of the final state photons are added to the electrons; iii) a truly 
multi-photon generator based on an exponentiated calculation (BHLUMI) [6,30]; it includes 
complete (9(a), O(aL) and 0(a 2 L 2 ) while 0(a 2 L) and 0(a 3 L 3 ) are incomplete; it generates 
explicitly 4-momenta of all photons above an arbitrary (user-defined) energy threshold, typi- 
cally a fraction k° (typically 10~ 4 ) of the beam energy. The Bhabha cross section calculated 
with BHLUMI will be compared to the one calculated with the hybrid calculation consisting 
of OLDBIS plus higher-order contributions from LUMLOG (LUMLOG#o)- The cross section 
differences BHLUMI - OLDBIS and BHLUMI - (OLDBIS + LUMLOG/™) are studied as 
a function of variations in the event selection parameters. Note that BHLUMI — OLDBIS is 
dominated by 0{a 2 L 2 ), G(a 2 L) and G(a 3 L 3 ) while BHLUMI - (OLDBIS + LUMLOG/™) is 
dominated by 0(a 2 L) and 0(a 3 L 3 ). 

Only the QED t-channel part of the generators is used, with photon vacuum-polarization 
switched off. We use an improved version of the BHLUMI event generator as discussed in 
Ref. [44]. BHLUMI — OLDBIS is used to estimate the higher-order contributions. We choose 
BHLUMI because the BHLUMI Monte Carlo distributions are in excellent agreement with 
the data distributions for all LEP experiments [47-52] A quantitative measurement of doubly 
radiative events [53] has shown consistency with the BHLUMI expectations and also with 
OLDBIS + LUMLOG#o expectations, while OLDBIS alone fails to describe this contribution, 
as expected. However, although the MC differential distributions agree with the data, the 
absolute scale of the integrated cross section remains uncertain, since the bulk of the radiative 
corrections are either virtual or involve soft (< 5 MeV) photons. 

In order to set the scale for the following numerical investigation let us remind the reader 
that the LEP1 experiments have reached in 1993-94 a systematic experimental uncertainty in 
the measuring the SABH luminosity cross section better than 0.10% [3-5]. 

2.6.1 Reference event selections 

We define an imaginary detector, consisting in a pair of cylindrical calorimeters covering the 
region between 62 and 142 mm radially out from the beam pipe centre and located at 2460 mm 
from the interaction point, at opposite sides of it. The beams are pointlike and centered within 
the beam pipe. The calorimeters are each divided into 32 azimuthal segments, subdivided into 
32 radial pads. A parton (electron or photon) deposits all its energy in the pad it hits. Photons 
and electrons from Bhabha events that hit the detector within a region of ±16 radial pads and 
5 azimuthal segments centered on the pad struck by the largest energy parton are combined 
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into a cluster. The cluster energy is the pad energy sum. Coordinates of the cluster centroid 
are the energy weighted average polar coordinates (i£, </>), summing over all pads in the cluster. 
Partons falling outside the principal cluster can originate secondary clusters, with no overlap. 
Only one cluster, the most energetic of all clusters, is used. Bhabha events are selected using 
the cluster energy E c i uster and the radial coordinate of the cluster centroid in both calorimeters. 

We then define a reference small-angle selection for Bhabha events (RSA selection). The 
radial acceptance edges for Bhabha events are set at pad boundaries. The "Wide" acceptance 
boundary extends up to two pads away from the detector inner and outer edges (27.236 < 
8 < 55.691 mrad). The "Narrow" acceptance boundary extends up to six pads away from 
the detector inner and outer edges (31.301 < 8 < 51.626 mrad). A similar angular range 
is covered by the OPAL, L3, ALEPH luminometers [54-56]. An event is selected when the 
cluster coordinates are within the Wide acceptance at one side (side 1) and within the narrow 
acceptance at the opposite side (side 2). Events must satisfy the criterion 0.5(a:i + x 2 ) > 0.75, 
with x = E c i uster / Eb eam . Selection criteria are also applied on the acoplanarity (0.2 rad) and 
the acollinearity (10 mrad) between the electron and positron clusters. 

Another selection is also considered, similar to the previous one but extending over the 
angular range covered by the DELPHI luminometer [57] (RLA: reference large-angle selection). 
The calorimeters are located at 2200 mm from the interaction point and cover radially the region 
between 6.5 and 41.7 cm. A cluster is formed starting from the most energetic particle hitting 
the calorimeter and considering all particles whose angular distance (A#, A</>) (in radians) from 
the initial one satisfies the two shower separation condition (determined from the comparison 
with the data) (A0/O.O3) 2 + (A</>/0.87) 2 < 1. The cluster energy is the sum over the energies of 
all particles inside the cluster, while the cluster coordinates are given by the energy weighted 
sum of their polar coordinates. Bhabha events are selected by cutting on the minimum cluster 
energy min(a;i, x 2 ), on the acoplanarity (20°) and on the cluster radial coordinate. The radial 
acceptance is defined on the Narrow side by the condition 43.502 < 8 < 113.151 mrad and on 
the Wide side by the condition 38.629 < 8 < 126.592 mrad. 

2.6.2 Comparison of exponentiated and order-by-order calculations 
First-Order Calculation 

The Bhabha cross-section for the RSA and RLA selections has been calculated with OLDBIS. 
The results are shown in figure 4 for the RSA selection, where the cross section is subdi- 
vided into a;-bins, separately for the narrow acceptance side and for the large acceptance side 
(^Narrow, %Wide)- A sample of 3 X 10 9 events is used. The total Bhabha cross section within the 
RSA acceptance is 75.589 ±0.009 nb. Displacing the generation minimum angle 8 9 ™ n from 10.4 
mrad as recommended in [6,30] to 5.2 mrad changes the accepted cross section by 0.0039(6) 
nbarns. No sizeable k° (=E^/ Eb eam ) dependence is observed when varying k° from 10~ 4 to 
10" 5 . 
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Figure 4: OLDBIS Bhabha cross section (nb) in 
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phase space bins for the RSA selection (see text). 



Higher-Order Leading-Log Contribution 

The cross section difference LUMLOG i/ o= LUMLOG(all orders) - LUMLOG(first order) is 
used to estimate the higher-order leading-logarithmic contribution (figure 5) for the RSA se- 
lection. In LUMLOG only the initial state radiation has an impact on the measured cluster 
energies and angles, because the 4-momenta of the final state photons are combined together 
with the electrons. A sample of 2.1 X 10 9 events is used. There is a total higher-order leading-log 
contribution of 0.144 ± 0.008 nb to the Bhabha cross section within the RSA acceptance: the 
higher-order contribution is negative in the phase-space region dominated by singly radiative 
events; it is positive in the non radiative Bhabha peak and in the phase-space region of hard 
doubly radiative events. 

Exponentiated Calculation 

The Bhabha cross-section in phase-space bins for the RSA selection obtained with BHLUMI 
is presented in figure 6. A sample of 1.6 X 10 9 events is used. The total Bhabha cross section 
accepted by the RSA selection is 75.712 ± 0.006 nb. The accepted cross section changes by 
< 10~ 5 when decreasing the t 9 ™ n (minimum generated four-momentum transfer squared) value 
as recommended in [6,30] to half of it. 

Comparison of Exponentiated and Order-by-Order Calculations 

The BHLUMI and OLDBIS cross sections differ for the RSA selection by (0.16±0.01)%, showing 
that the estimated contribution to the accepted cross section from higher-order radiative effects 
is very small. This estimate is also in reasonable agreement with the LUMLOG#o expectation 
of (0.19 ±0.01)%. 

A similar study for the RLA selection results in a BHLUMI — OLDBIS relative difference 
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Figure 5: LUMLOG higher-order contribution to the Bhabha cross section (nb) in phase space bins for 
the RSA selection (see text). 
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Figure 6: BHLUMI Bhabha cross section (nb) in phase space bins for the RSA selection (see text). 



of (-0.08 ± 0.01)% to be compared with a LUML0G# o expectation of (-0.03 ± 0.01)%. 



18 



2.6.3 Dependence on energy and acollinearity cuts 



The cross section relative difference (BHLUMI - OLDBIS)/BHLUMI(RSA), where BHLUMI(RSA) 
refers to the RSA selection, is studied in table 8 for several selection criteria on energy and 
acollinearity. With x™j n we mean that the energy cut min(xi, x 2 ) > x™j n 1S applied. Through 
transverse momentum conservation, energy and acollinearity cuts are strongly correlated in 
events with initial state radiation. The relative difference BHLUMI — OLDBIS is indicative 
of the higher-order contribution, which clearly appears in table 8 to be huge for large x^ n . It 
becomes progressively smaller for smaller x^ n . It should be stressed that the h.o. corrections 
are small (at the per mille level) over a very broad region of x™j n and acollinearity. 



A second estimate of the Bhabha cross section with higher-order radiative corrections can 
be obtained with OLDBIS + LUML0G# O - The three generator relative difference (BHLUMI 
- (OLDBIS + LUML0G# o ))/BHLUMI(RSA) in table 8 shows that the h.o. corrections in 
BHLUMI and in LUMLOG track each other very well, giving confidence that the h.o. contri- 
butions are in fact small when they are estimated to be so. The unstable region is limited to 
very large x%$ n . The BHLUMI and OLDBIS + LUML0G# o Bhabha cross sections agree at 
the 0.1% level over an extremely broad range of energy and acollinearity cuts. 

The cross section differences BHLUMI - OLDBIS and BHLUMI - (OLDBIS + LUMLOG/™) 
for the RSA selection change by (—0.013 ± 0.009)% when the acoplanarity cut is not applied. 

For the RLA selection the cross section differences BHLUMI - OLDBIS and BHLUMI - 
(OLDBIS + LUMLOG#o) normalized to the BHLUMI result are shown in table 9 as a function 
of the cut on x^ n . The higher-order contribution to the Bhabha cross section for the RLA 
selection both in BHLUMI and in LUMLOG is very small over a broad range of x c r 



„cut 



BHLUMI— OLDBIS 



BHLUMI-(OLDBIS+LUMLOG#o) 





Acollinearity cut (rad) 


T cut 
min 


0.005 0.010 no cut 


0.999 


11.35(1)% 10.86(1)% 10.61(1)% 


0.99 


4.65(1)% 4.45(1)% 4.35(1)% 


0.90 


0.69(2)% 0.60(1)% 0.58(1)% 


0.85 


0.68(1)% 0.25(1)% 0.24(1)% 


0.75 


0.75(1)% 0.12(1)% -0.00(1)% 


triang. 


0.78(1)% 0.16(1)% -0.09(1)% 


0.50 


0.83(1)% 0.26(1)% 0.06(1)% 





Acollinearity cut (rad) 


T cut 
min 


0.005 0.010 no cut 


0.999 


2.19(2)% 2.10(1)% 2.05(1)% 


0.99 


0.98(2)% 0.94(2)% 0.92(2)% 


0.90 


0.19(2)% 0.15(2)% 0.14(2)% 


0.85 


0.15(2)% 0.06(2)% 0.06(2)% 


0.75 


0.13(2)% -0.00(2)% -0.03(2)% 


triang. 


0.12(2)% -0.03(2)% -0.09(2)% 


0.50 


0.18(2)% -0.02(2)% -0.07(2)% 



Table 8: Cross section differences BHLUMI-OLDBIS and BHLUMI-(OLDBIS+LUMLOG i/ o) normalized 
to the BHLUMI Bhabha cross section for the RSA selection. The label "triangular" stands for the cut 
0.5(a:i + x 2 ) > 0.75. 
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cut 
min 


BHL-OB 


0.9 


0.76(1)% 


0.8 


0.l0(l)% 


0.7 


-0.06(l)% 


0.6 


-0.08(l)% 


0.5 


-0.05(l)% 



cut 
min 


BHL— (0B+LL#o) 


0.9 


0.2l(l)% 


0.8 


0.03(l)% 


0.7 


-0.03(l)% 


0.6 


-0.05(l)% 


0.5 


-0.06(l)% 



Table 9: Cross section differences BHLUMI-OLDBIS and BHLUMI-(0LDBIS+LUML0G# o ) normalized 
to the BHLUMI Bhabha cross section for the RLA selection. 



2.6.4 Wide- Wide, Narrow-Narrow versus Wide-Narrow acceptance 



In the reference selections (RFA and RLA) an asymmetric acceptance (Wide on one side and 
Narrow on the opposite side) is used. All 4 LEP experiments use an asymmetric acceptance for 
the LEP luminosity measurement. We study in table 10 how the results change when using a 
symmetric (Wide- Wide or Narrow-Narrow). The BHLUMI — OLDBIS cross section difference 
becomes large (0.77(l)% for the Narrow-Narrow acceptance). A similar result is also obtained 
using LUMLOGffo and then the BHLUMI - (OLDBIS + LUMLOG^o) difference is small. 
We thus conclude that the higher-order contributions to the accepted Bhabha cross section, as 
estimated with BHLUMI or LUMLOG, are largely reduced when using an asymmetric Wide- 
Narrow acceptance. 





WN WW NN 


BHLUMI 


75.7l2(5)n6 H7.9l8(6)n6 73.344(5)n6 


OLDBIS 


75.589(8)n6 H7.2l9(9)n6 72.78l(8)n6 


LUMLOGffo 


0.l44(8)n6 0.568(9)n6 0.465(8)n6 


(BHL— OB)/BHL 


0.l6(l)% 0.59(l)% 0.77(l)% 


(BHL— OB— LL# )/BHL 


-0.03(2)% 0.ll(2)% 0.13(2)% 



Table 10: Comparison of BHLUMI, OLDBIS and LUML0G# o Bhabha cross sections for Wide-Narrow, 
Wide-Wide, Narrow-Narrow event selections. All other cuts as in the RSA selection. 



2.6.5 Multiple photon radiation 



A very relevant property of exclusive exponentiation is that there are many more multi-photon 
events than expected from perturbation theory at a fixed order in a. In a sample of 10 6 BHLUMI 
Bhabha events, the events have up to eight photons with energy larger than k°Eb eam (tt 5 MeV), 
as shown in figure 7. This may enhance the difference between cross section calculations 
performed with BHLUMI and with OLDBIS + LUML0G# O - In the following we study the 
stability of the BHLUMI - OLDBIS and BHLUMI - OLDBIS - LUMLOG(ho) differences in 
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table 8 and in table 9 when varying those parameters in the experimental selection which are 
sensitive to the presence of many photons. 

Lower Energy Photon Cut-off 

We define a K c parameter (in MeV) expressing the sensitivity to soft photons: the detector 
is fully efficient for photons of energy larger than K c . An implicit K c cut-off is present in 
BHLUMI at K c = k°Eb eam (5 MeV) for the cross sections calculations presented above. The 
relative variation of the BHLUMI Bhabha cross section when varying K c is reported in table 
11 for the RSA selection and in table 12 for the RLA selection. The effect is at most of 
—0.03% for the RSA acceptance in the extreme case of K c =500 MeV. The relative changes 
in the BHLUMI and OLDBIS cross sections are compared in figures 8 and 9. The large-x 
region is very different; most of the difference has already disappeared for x^ n =0.9. LUMLOG 
remains unaffected: it has in the output only the electron and positron 4-momenta with the 
final state photons combined with the electrons/positrons. Hence, the effect on the relative 
cross section differences BHLUMI - OLDBIS and BHLUMI - (OLDBIS + LUMLOG/™) for 
the RSA selection is at most -0.030(4)%. 
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Figure 7: Distribution in number of emitted photons for a sample of 10 6 unweighted BHLUMI events. 
(The Removal flag is switched on in BHLUMI, with K c = k°Eb eam = 5 MeV). 

Cluster Size 

The relative variation of the accepted Bhabha cross section with respect to the RSA selection 
when changing the cluster size is studied in figure 10 using BHLUMI generated events and using 
OLDBIS generated events. For large cluster sizes BHLUMI and OLDBIS track each other very 
well and the BHLUMI — OLDBIS relative difference observed for the RSA selection remains 
unchanged. On the contrary, for small cluster sizes, the effect of many photons in BHLUMI 
generated events shows up strongly. The LUMLOG result remains unaffected. Thus, for the 
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cut 
min 


K c (MeV) 


1U OU 1UU ouu 


0.999 


-0.025(4)% -0.75(2)% -3.51(5)% -9.45(8)% 


0.90 


-0.001(1)% -0.002(1)% -0.004(2)% -0.029(4)% 


0.85 


< 10" 5 -0.003(1)% -0.003(1)% -0.012(3)% 


triangular 


-0.004(2)% -0.015(3)% -0.018(3)% -0.030(4)% 



Table 11: Variation of the BHLUMI Bhabha cross section when changing the photon minimum detectable 
energy K c . Normalization is with respect to the RSA selection with K c = k°E}, eam = 5 MeV. The label 
"triangular" stands for the cut 0.5(a;i + x^) > 0.75. 



cut 
min 


K c (MeV) 


10 50 100 500 


0.9 


-0.0005(2)% -0.0032(5)% -0.0067(7)% -0.033(2)% 


0.7 


-0.0003(2)% -0.0010(3)% -0.0015(3)% -0.0085(8)% 


0.5 


< 10" 6 -0.0002(1)% -0.0005(2)% -0.0025(5)% 



Table 12: Variation of the BHLUMI Bhabha cross < 
energy K c from K c = 5 MeV for the RLA selection 




0.50 0.65 0.75 0.85 1 



ion when changing the photon minimum detectable 



-0.072(7) 


-0.111(8) 


3.79(5) 


-9.45(8) 


0.072(7) 


0.50(2) 


1.11(3) 


3.65(5) 


0.009(2) 


0.080(7) 


0.54(1) 


-0.137(9) 


< 0.001 


0.010(3) 


0.105(8) 


-0.106(8) 



0.85 0.90 0.99 0.999 1 



Figure 8: Percentage variation of the BHLUMI Bhabha cross section when setting the photon minimum 
detectable energy K c to 500 MeV (see also figure 3) instead of K c = 5 MeV. 
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< 0.001 
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< 0.001 
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0.99 



0.999 



Figure 9: Percentage variation of the OLDBIS Bhabha cross section when setting the photon minimum 
detectable energy K c to 500 MeV (see also figure 1) instead of K c = 5 MeV. 



RSA selection we can exclude an effect larger than ±0.007% on the BHLUMI — OLDBIS and 
on the BHLUMI - (OLDBIS + LUML0G# o ) cross section differences. 

Cluster Coordinate 

The energy weighting algorithm for extracting the cluster coordinates couples the coordinates 
to the cluster size. A different coordinate reconstruction algorithm (PADMAX) is then used: 
we select the pad with the largest energy deposit and use the 4-momentum sum of the partons 
which enter that pad to calculate an impact point in the pad; the impact point so calculated 
defines the cluster coordinates, independent of the cluster dimensions. The BHLUMI cross 
section when changing from 6 (</>) energy weighted coordinates to PADMAX coordinates in the 
RSA selection changes by (—0.088 ± 0.003)%. The OLDBIS cross section when changing from 8 
(</>) energy weighted coordinates to PADMAX coordinates changes by (—0.091 ±0.005)%. The 
LUMLOG result is unaffected. The effect on the BHLUMI - OLDBIS and on the BHLUMI 
— (OLDBIS ± LUMLOG#o) cross section differences in the RSA selection when using the 
PADMAX coordinates instead of the energy weighted coordinates is (0.003 ± 0.006)%. 



2.6.6 Summary 

We have shown that there is a strong correlation between the magnitude of the 0(a 2 ) radiative 
corrections to the Bhabha cross section and distinctive characteristics of the experimental 
Bhabha event selection. In particular, we have shown that the Bhabha selections used by 
the LEP experiments to measure the accelerator luminosity minimize the sensitivity to 0(a 2 ) 
radiative corrections. 
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Figure 10: Relative variation of the accepted Bhabha cross section with respect to the RSA selection when 
changing cluster radial (PAD's) and azimuthal (SEGments) dimensions. A cluster extends for ±Apad pads 
and ±Nseg segments around the pad containing the largest energy deposit. A pad subtends a polar angle of 
about 1 mrad; a segment covers azimuthally an angle of 11.25 degrees. The RSA selection has Apad = 16 
and N S eg = 2. 

The 0(a 2 ) contributions have been estimated using BHLUMI - OLDBIS and LUMLOG#o = 
LUMLOG a M_ or< iers — LUMLOG/; rS 4_ or( ier- The cross section differences BHLUMI — OLDBIS 
and BHLUMI - (OLDBIS + LUML0G# o ) are very small (at the per mille level) in a broad 
region of phase space around the experimental selections. We have considered two angular 
ranges 27 < 8 < 57 mrad and 44 < 8 < 113 mrad, with a variety of energy and acollinearity 
cuts. The sensitivity to the possible presence of many photons, predicted by exclusive expo- 
nentiation, the effect of small or large cluster sizes and different ways of reconstructing the 
cluster coordinates have been investigated. Large cluster sizes, rather soft energy cuts and a 
Wide-Narrow method are very effective in minimizing the cross section differences BHLUMI — 
OLDBIS and BHLUMI - (OLDBIS + LUML0G# o ). Vice versa, these same effects could be 
used to enhance the sensitivity to the 0(a 2 ) radiative corrections in order to perform measure- 
ments and test the theory predictions. 
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2.7 Comparisons of event generators for small-angle Bhabha scat- 
tering 



In contrast to the previous section, where we have seen results from many variants of ES's 
with varying cut parameters but for only three types of QED calculations, here we shall limit 
ourselves to "only" four ES's (two of which very close to realistic experimental situations), 
but we shall discuss all the available theoretical calculations. The outline of this section is 
the following: the actual comparisons will be presented first at the 0(a 1 ) level, in order to 
determine the basic technical precision, and later for more advanced QED matrix elements 
beyond 0(a 1 ), in order to explore physical precision. These comparisons will be done first for 
LEP1 energy and later will be also extended to LEP2 energies. 



Forward hemisphere 



BARE1 



Backward hemisphere 




WW:^(C.C). NN:0i<E(ft 



3iV 



NW: 9 1 G (Cm, Cax), 62 e (Cn,CJ 



QiV 



mm' maxi 



S > U r 



Figure 11: Geometry and acceptance of the simple (non-calorimetric) ES BARE1. This ES restricts polar 
angles 9{ in the forward/backward hemispheres and requires a certain minimum energy to be detected 
simultaneously in both hemispheres. Photon momentum is not constrained at all. The entire "fiducial" 



0-range, i.e. wide (W) range, is (Cm, Cax) 



(0.024, 0.058) rad and the narrow (N) range is (6^ in , 6>£ ax ), 



where Cn = Cin + Cx = Cax " Se and S = (Cax " Cin)/16. This ES can be symmetric Wide- 
Wide (WW) or Narrow-Narrow (NN), or asymmetric Narrow-Wide (NW), see the description in the figure. 
The energy cut s' > u m i n s involves momenta of outgoing e^ (s' = (q + + Q~) 2 ) only. 



2.7.1 Event selections 

One cannot talk about the cross section for the small-angle Bhabha (SABH) process without 
defining precisely all cuts, or, in other terms, without specifying the ES. The most interesting 
ES is that of the actual experiment. LEP1 and LEP2 experiments employ in the measurement 
of the small-angle Bhabha scattering cross section a rich family of ES's. They do, however, 
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Forward hemisphere 



CAL01 



Backward hemisphere 




fid 2x 
WW: Of e (Cin,Cax), NN: 9f e (Cin,Cx) 

NW:^e(Cin ? eax), ^e(Cin,Cx) 



T?cl T?d \ - 7Tl2 

-C/i -C/9 > Z m i n Jl,, 



Figure 12: Geometry and acceptance of the calorimetric ES CAL01. This ES restricts polar angles 9{ in 
the forward/backward hemispheres and requires a certain minimum energy to be detected simultaneously in 
both hemispheres. The entire "fiducial" #-range, i.e. wide (W) range, is (^ n ,Cax) = (0.024,0.058) rad 
and the narrow (N) range is (0^,_CJ, where 9% in = Cin + h, Cx = Cax " h and ^ = (Cax " 
#^i n )/16. This ES can be symmetric Wide-Wide (WW) or Narrow-Narrow (NN), or asymmetric Narrow- 
Wide (NW), see the description in the figure. The energy cut involves the definition of the cluster, the 
cluster center (9f,<j)f), i = 1,2, is identical to the angular position of the positron in the forward and 
the electron in the backward hemisphere. The angular "cone" of radius 6 = 0.010 rad around e^ is called 
cluster. The cone/cluster in the 9, (f> plane is an elongated ellipsis, due to smallness of theta. The total 
energy registered in the cluster is denoted by Ef. (Note that (j>i = for back-to-back configuration.) 



have essential common features. The most important is the "double tag". It means that e + 
and e~ are both detected with a certain minimum energy and minimum scattering angle in the 
forward and backward direction, close to the beams. The other important feature of the typical 
experimental ES is that (except for rare cases) the photons and e ± cannot be distinguished - 
only the combined energy and angle is registered. It is said that the typical experimental ES 
is calorimetric. On the other hand, for comparing theoretical calculations it is useful to deal 
with simplified ES's, in which only e ± are measured and the accompanying bremsstrahlung 
photons (e ± pairs) are ignored. The "double tag" is done on "bare e ± " . Actually, in order 
to compare efficiently numerical results from the various programs, we employed the family 
of four ES's connecting in an almost continuous way the experimentally unrealistic (but use- 
ful for theorists) examples of ES's to experimentally realistic (but difficult for some class of 
theoretical calculations) ones. In order to compare theoretical results for SABH, we use one 
simple non-calorimetric ES called BARE1, see Figs. (11), and three calorimetric ES's called 
CALOl, CAL02 and SICAL2, with increasing degrees of sophistication. They are defined in 
Figs. (12,13) and Fig. (14). The last one, SICAL2 of Fig. (14), corresponds very closely to the 
ES of the real silicon detector of OPAL or ALEPH. 
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Forward hemisphere 



CAL02 




Backward hemisphere 



<t>2 



<J>? 2x 

WW: Of e (Cin,Cax), NN: Of G (CCx) 
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Figure 13: Geometry and acceptance of the calorimetric ES CAL02. This ES restricts polar angles 9{ in 
the forward/backward hemispheres and requires a certain minimum energy to be detected simultaneously 
in both hemispheres. The entire "fiducial" #-range, (^ n) ^ax) = (0.024,0.058) rad, includes the wide 
(W) range (C in , Cax) and the narrow (N) range (C n , Cx), where Cin = *Ln + Cax = tx " S e , 
Se = - ^ in )/16, and Cn = + ™e, Cx = tx " ^e- This ES can be symmetric Wide-Wide 

(WW) or Narrow-Narrow (NN), or asymmetric Narrow-Wide (NW), see the description in the figure. The 
energy cut involves the definition of the cluster: the cluster center [0f,<j>f), i = 1,2, is identical to the 
angular position of the positron in the forward and electron in the backward hemisphere. The angular 
"plaquette" (6f + l.58 e , 6f - l.58 e ) X {<f>f + 1.5fy, <f>f - 1.5fy), where ^ = 2x/32, around e ± is called 
cluster. The total energy registered in the cluster is denoted by Ef. (Note that (f>i = (f>2 for back-to-back 
configuration.) 



2.7.2 First order - technical precision 

We start the numerical comparisons of the various theoretical calculations with the calibration 
exercise in which we limit ourselves to strict 0(a 1 ) with Z exchange, up-down interference 
and vacuum polarization switched off, i.e. we examine pure photonic corrections without up- 
down interferences. We calculate the corresponding total cross section for all our four ES's 
at the LEP1 energy, y^s = 92.3 GeV. The purpose of this exercise is to eliminate possible 
trivial normalization problems in the core MC programs and in the testing programs which 
implement our ES's. Since 0(a 1 ) is unique and common, the difference of the results will 
be entirely due to numerical/technical problems and, following ref. [11] where the analogous 
exercise of this type was done for the first time, we call it the "technical precision" of the 
involved calculations/programs. The results are shown in Tab. 13. Since tables are hard to 
read, we always include a figure which contains exactly the same result in the pictorial way. 
In the figure, one of the cross sections is used as a reference cross section and is subtracted 
from the other ones. It is plotted however on the horizontal line with its true statistical error. 
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Figure 14: Geometry and acceptance of the calorimetric ES SICAL2. This ES restricts polar angles 9{ in 
the forward/backward hemispheres and requires a certain minimum energy to be detected simultaneously 
in both hemispheres. No restrictions on azimuthal angles <f>i are there. The entire "fiducial" #-range, 

(^min^Lx) = (0.024,0.058) rad, includes the wide (W) range (Cin>Cax) and the narrow (N) range 
(^min' ^max) exactly as depicted in the figure. This ES can be symmetric Wide-Wide (WW) or Narrow- 
Narrow (NN), or asymmetric Narrow-Wide (NW). The energy cut and #-cuts involve the definition of the 
cluster. Eeach side detector consists of 16 X 32 equal plaquetes. A single plaquete registers the total energy 
of electrons and photons. The plaquete with the maximum energy, together with its 3 X 3 neighborhood, is 
called cluster. The total energy registered in the cluster is Ef and its angular position is [0f,<j>f), i = 1,2. 
More precisely the angular position of a cluster is the average position of the centers of all 3 X 3 plaquetes, 
weighted by their energies (the definitions of (f>'s are adjusted in such a way that (f>i = fa for back-to-back 
configuration). The plaquetes of the cluster which spill over the angular range (outside thick lines) are also 
used to determine the total energy and the average position of the cluster (see backward hemisphere). 



Here Tab. 13 is visualized in Fig. 15. In this figure, the cross sections from the Monte Carlo 
OLDBIS (an improved version of the MC program written originally by Berends and Kleiss in 
PETRA times, now part of BHLUMI) is used as a reference. As we see, all calculations agree 
well within 3 X 10~ 4 relative deviation. The apparent discrepancy of the 0(a 1 ) SABSPV for 
the SICAL2 ES is not statistically significant. The cross section from the non-Monte-Carlo 
type of calculation NLLBHA is available only for the simplest BARE1. As we have already 
discussed, the photonic radiative corrections for the SABH process scale smoothly with energy, 
so we regard this test to be valid for LEP2 energies within a factor two, i.e. within 6 X 10~ 4 . 



2.7.3 Beyond first order - physical precision 



Having found good agreement of the various calculations at the first order level, we now reinstall 
the photonic corrections beyond first order. More precisely we keep again Z exchange, up-down 
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•^min 


OLDBIS [nb] 


SABSPV [nb] 


BHAGEN95 [nb] 


NNLBHA [nb] 


BHLUMI [nb] 




(a) BARE1 


.100 


166.079 ± .013 


166.070 ± .024 


.000 ± .000 


166.070 ± .017 


166.046 ± .021 


.300 


164.772 ± .013 


164.762 ± .012 


164.756 ± .012 


164.767 ± .016 


164.740 ± .021 


.500 


162.277 ± .013 


162.263 ± .012 


162.258 ± .012 


162.265 ± .016 


162.241 ± .021 


.700 


155.465 ± .013 


155.452 ± .012 


155.444 ± .012 


155.453 ± .015 


155.431 ± .020 


.900 


134.417 ± .012 


134.401 ± .023 


134.394 ± .012 


134.393 ± .014 


134.390 ± .020 




(b) CAL01 


.100 


166.361 ± .013 


166.353 ± .024 


.000 ± .000 


.000 ± .000 


166.329 ± .021 


.300 


166.081 ± .013 


166.071 ± .021 


166.074 ± .013 


.000 ± .000 


166.049 ± .021 


.500 


165.319 ± .013 


165.311 ± .012 


165.312 ± .013 


.000 ± .000 


165.287 ± .021 


.700 


161.823 ± .013 


161.817 ± .024 


161.818 ± .013 


.000 ± .000 


161.794 ± .021 


.900 


149.942 ± .013 


149.934 ± .023 


149.934 ± .013 


.000 ± .000 


149.925 ± .020 




(c) CAL02 


.100 


131.061 ± .012 


131.070 ± .022 


131.051 ± .010 


.000 ± .000 


131.032 ± .019 


.300 


130.769 ± .012 


130.778 ± .022 


130.758 ± .010 


.000 ± .000 


130.739 ± .019 


.500 


130.206 ± .012 


130.214 ± .022 


130.194 ± .010 


.000 ± .000 


130.176 ± .019 


.700 


127.555 ± .012 


127.565 ± .022 


127.546 ± .010 


.000 ± .000 


127.528 ± .019 


.900 


117.557 ± .011 


117.572 ± .025 


117.543 ± .010 


.000 ± .000 


117.541 ± .018 




(d) SICAL2 


.100 


132.011 ± .012 


131.965 ± .023 


132.004 ± .028 


.000 ± .000 


131.984 ± .019 


.300 


131.900 ± .012 


131.862 ± .021 


131.893 ± .027 


.000 ± .000 


131.872 ± .019 


.500 


131.587 ± .012 


131.539 ± .018 


131.581 ± .027 


.000 ± .000 


131.559 ± .019 


.700 


128.363 ± .012 


128.306 ± .016 


128.364 ± .027 


.000 ± .000 


128.338 ± .019 


.900 


117.843 ± .011 


117.795 ± .012 


117.811 ± .027 


.000 ± .000 


117.828 ± .018 



Table 13: Monte Carlo results for the symmetric Wide-Wide ES's BARE1, CAL01, CAL02 and SICAL2, 
for the 0(a 1 ) matrix element. Z exchange, up-down interference and vacuum polarization are switched 
off. The center of mass energy is ^/s = 92.3 GeV. Not available x-sections are set to zero. 



interference and vacuum polarization switched off, but compare numerical results which include 
0(a 2 L 2 ), 0(a 2 L) and 0(a 3 L 3 ) contributions due to photon bremsstrahlung. We do not include 
production of light fermion pairs unless stated otherwise. The numerical results are shown in 
Tab. 14 and Fig. 16. In the figure, the cross section from the second order exponentiated 
Monte Carlo BHLUMI is used as a reference cross section. The differences between various 
calculations now represent not only technical precision, but also physical precision because the 
cross sections are calculated using different QED matrix elements. 

The results shown in Tab. 14 and Fig. 16 have remarkable properties. For values of the 
energy-cut variable in the experimentally interesting range 0.25 < z m i n < 0.75, the cross sec- 
tion from the programs BHLUMI and SABSPV agree throughout all the four ES's, from the 
unrealistic BARE1 to very realistic SICAL2, to within 1.0 X 10~ 3 relative deviation. This agree- 
ment is definitely better than the difference between BHLUMI and OLDBIS+LUMLOG, which 
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Figure 15: Monte Carlo results for the symmetric Wide-Wide ES's BARE1, CALOl, CAL02 and SICAL2, 
for the 0(a 1 ) matrix element. Z exchange, up-down interference and vacuum polarization are switched 
off. The center of mass energy is ^ = 92.3 GeV. In the plot, the cross section from the program OLDBIS 
(part from BHLUMI 4. 02. a, originally written by Berends and Kleiss) is used as a reference cross section. 



in the last years was routinely used (see Refs. [6,58]) in order to estimate missing higher order 
and subleading corrections. Remarkably, the OLDBIS+LUMLOG results coincide extremely 
well with BHAGEN95. Let us note that the OLDBIS+LUMLOG matrix element does not ex- 
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BHLUMI [nb] 


SABSPV [nb] 


BHAGEN95 [nb] 


OBI+LMG [nb] 


NLLBHA [nb] 




(a) BARE1 


.100 


166.892 ± .006 


166.795 ± .028 


.000 ± .000 


166.672 ± .017 


166.948 ± .000 


.300 


165.374 ± .006 


165.323 ± .028 


165.190 ± .012 


165.187 ± .017 


165.448 ± .000 


.500 


162.530 ± .006 


162.529 ± .028 


162.330 ± .012 


162.365 ± .017 


162.581 ± .000 


.700 


155.668 ± .006 


155.751 ± .026 


155.466 ± .012 


155.519 ± .017 


155.617 ± .000 


.900 


137.342 ± .006 


137.528 ± .022 


137.188 ± .011 


137.210 ± .017 


137.201 ± .000 




(b) CALOl 


.100 


167.203 ± .006 


167.106 ± .028 


.000 ± .000 


167.000 ± .017 


.000 ± .000 


.300 


166.795 ± .006 


166.715 ± .028 


166.618 ± .012 


166.623 ± .017 


.000 ± .000 


.500 


165.830 ± .006 


165.768 ± .014 


165.661 ± .014 


165.686 ± .017 


.000 ± .000 


.700 


162.237 ± .006 


162.203 ± .027 


162.048 ± .014 


162.053 ± .017 


.000 ± .000 


.900 


151.270 ± .006 


151.272 ± .025 


150.823 ± .014 


150.707 ± .017 


.000 ± .000 




(c) CAL02 


.100 


131.835 ± .006 


131.755 ± .027 


131.658 ± .007 


131.632 ± .016 


.000 ± .000 


.300 


131.450 ± .006 


131.393 ± .027 


131.285 ± .012 


131.274 ± .016 


.000 ± .000 


.500 


130.727 ± .006 


130.708 ± .027 


130.575 ± .012 


130.584 ± .016 


.000 ± .000 


.700 


127.969 ± .006 


127.999 ± .027 


127.802 ± .014 


127.802 ± .016 


.000 ± .000 


.900 


118.792 ± .006 


118.879 ± .029 


118.293 ± .013 


118.201 ± .015 


.000 ± .000 




(d) SICAL2 


.100 


132.816 ± .006 


132.612 ± .026 


132.611 ± .028 


132.582 ± .016 


.000 ± .000 


.300 


132.553 ± .006 


132.427 ± .025 


132.420 ± .028 


132.405 ± .016 


.000 ± .000 


.500 


131.985 ± .006 


131.966 ± .022 


131.962 ± .027 


131.965 ± .016 


.000 ± .000 


.700 


128.672 ± .006 


128.691 ± .019 


128.620 ± .027 


128.610 ± .016 


.000 ± .000 


.900 


119.013 ± .006 


119.075 ±.015 


118.561 ± .027 


118.488 ± .015 


.000 ± .000 



Table 14: Monte Carlo results for the symmetric Wide-Wide ES's BARE1, CALOl, CAL02 and SICAL2, 
for matrix elements beyond first order. Z exchange, up-down interference and vacuum polarization are 
switched off. The center of mass energy is ^/s = 92.3 GeV. Not available x-sections are set to zero. 



ponentiate properly 0(a 2 L) corrections, i.e. they are wrong in the soft photon limit. This may 
explain why BHLUMI and SABSPV, which do not have such problems, agree better. According 
to the authors, BHAGEN95 does not suffer of the same problem as it has the soft photon limit 
properly treated by construction, but some corrections are expected due to the approximate 
treatment of two hard photon emission. The result from NLLBHA is present only for unreal- 
istic BARE1 selection, and for 0.25 < z min < 0.75 it agrees to within 0.1% with BHLUMI and 
SABSPV. It is an interesting result because NLLBHA features complete 0(a 2 L) corrections, 
while all the other programs have only incomplete 0(a 2 L) contributions. In Tab. 14 and Fig. 16 
the results of BHLUMI, SABSPV and BHAGEN95 include exponentiation, and therefore they 
include necessarily 0(a 3 L 3 ) effects (incomplete). We therefore compare them with a version 
of NLLBHA which includes, besides 0(a 2 L), also 0(a 3 L 3 ) corrections. All the above results 
will be used as an input in our final estimate of the total theoretical uncertainty of SABH cross 
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Figure 16: Monte Carlo results for the symmetric Wide-Wide ES's BARE1, CAL01, CAL02 and SICAL2, 
for matrix elements beyond first order. Z exchange, up-down interference and vacuum polarization are 
switched off. The center of mass energy is = 92.3 GeV. In the plot, the 0{a 2 )^ p s cross section c BHL 
from BHLUMI 4. 02. a is used as a reference cross section. 

section for LEP1/LEP2 energies. 

Finally, we present similar numerical comparisons of the calculations beyond 0(a 1 ) at one 
LEP2 energy y/s = 176 GeV. As before, since the tables are hard to read, we accompany 
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BHLUMI [nb] 


SABSPV [nb] 


BHAGEN95 [nb] 


OBI±LUM [nb] 




(a) CAL02 LEP2 


.100 


36.123 ±.003 


36.096 ± .008 


.000 ± .000 


36.060 ± .006 


.300 


36.013 ±.003 


35.992 ± .008 


35.963 ± .005 


35.958 ± .006 


.500 


35.807 ± .003 


35.796 ± .008 


35.762 ± .005 


35.761 ± .006 


.700 


35.001 ± .003 


35.005 ± .008 


34.951 ± .005 


34.948 ± .006 


.900 


32.324 ± .003 


32.341 ± .008 


32.173 ± .006 


32.145 ± .006 




(b) SICAL2 LEP2 


.100 


36.394 ± .003 


36.337 ±.011 


.000 ± .000 


36.322 ± .006 


.300 


36.316 ±.003 


36.284 ± .010 


36.271 ± .009 


36.270 ± .006 


.500 


36.150 ±.003 


36.147 ± .009 


36.139 ± .009 


36.142 ± .006 


.700 


35.193 ±.003 


35.203 ± .008 


35.173 ± .009 


35.171 ± .006 


.900 


32.383 ± .003 


32.405 ± .006 


32.243 ± .009 


32.224 ± .006 
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Table 15: In this table/figure we show cross sections for LEP2 center of mass energy, ^ = 176 GeV. 
Monte Carlo results are shown for various symmetric Wide-Wide ES's and matrix elements beyond first 
order. Z exchange, up-down interference and vacuum polarization are switched off. Not available x-sections 
are set to zero. In the plot, the 0(a 2 ) 
cross section. 



exp cross section c BHL from BHLUMI 4. 02. a is used as a reference 



the table with a figure which shows the same numerical result in a pictorial way (the caption 
is common for the table and figure). This way of presenting results in the form of the twin 
table/figure will be used often in the following. As before, in the figure one of the cross sections 
is used as a reference cross section and is subtracted from the other ones. The main result is 
shown in table/figure 15. Here, results are shown for the symmetric Wide- Wide variant of the 
CAL02 and SICAL2 ES's. As expected, the difference between the programs is almost the 
same! The higher order corrections at LEP2 are only slightly stronger. This result was already 
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anticipated when analyzing "scaling rules" derived from Tab. 2. From the scaling rules we also 
know that this result will be essentially the same for the wider angular range 3° — 6°. The 
practical message is that, within 20-30%, the precision estimates derived from the numerical 
exercises for the SABH process at LEP1 should be valid also for LEP2. 

Precision requirements at LEP2 are less stringent. In the figure, we draw a LEP2-type box 
which spans over 0.2% and extends over the experimentally interesting range 0.25 < z m i n < 
0.75. All programs come together within the above range. The above 0.2% limit will be used 
as an input in our final estimate of the total theoretical uncertainty of the SABH cross section 
for LEP2 energies. This limit has obviously a large safety margin, close to a factor of two. 

2.7.4 Asymmetric and very narrow event selections 

The numerical comparisons shown in the previous section were done, for pure technical reasons 
(less chances for programming errors in the testing programs), for the symmetric Wide- Wide 
version of the ES. As we know very well (see the introduction), the higher order contribu- 
tions are sensitive to the "asymmetricity" of the ES. In order to avoid any danger due to the 
above simplification, we have done another series of comparisons of the various calculations 
for the symmetric Narrow-Narrow and asymmetric Narrow- Wide versions of the ES's CAL02, 
which are defined in Fig. 13. Let us remind the reader that the variation of the difference 
BHLUMI-(OLDBIS+LUMLOG i/ o) over the WW, NN and NW selection was the cornerstone 
of the previous estimates [6,58] of the size of uncontrolled higher order photonic corrections (to- 
gether with technical precision). We believe that CAL02 is close enough to our most realistic 
ES SICAL2 and the results obtained for CAL02 are valid for SICAL2. Let us also recall that 
the typical experimental ES is of the asymmetric Narrow- Wide type. The corresponding results 
are shown in table/figure 16 for the matrix elements in the 0(a 2 ) class (we have checked that 
for the 0(a 1 ) level the same programs agree better than 0.03%, but we omit the corresponding 
table/plot due to lack of space). 

As we see in tables /figures 16 and 14, for all the three types of the CAL02 ES (WW, NN 
and NW), BHLUMI and SABSPV stay within 0.1% from one another for all the values of the 
energy-cut variable in the experimentally interesting range 0.25 < z m i n < 0.75. This is a new 
nontrivial result, which will be exploited to decrease the estimated error due to the higher order 
photonic corrections from 0.15% down to 0.1%. In a sense, we replace the old estimate based on 
BHLUMI - (OLDBIS + LUMLOG^o) with a new one based on BHLUMI-SABSPV. Hybrid 
Monte Carlo's (OLDBIS + LUMLOG^o) and BHAGEN95 are off of about 0.2% in the NN 
case but, noticeably, they are on the same ground as BHLUMI and SABSPV for the most 
interesting NW case. The above exercise was done for the LEP1 energy, and in view of the 
results shown in table/figure 15 and our "scaling rules" (see the introduction), we do not foresee 
any problem with extending its validity to LEP2 energies. 

As we already stressed in the introduction, for the purpose of LEP2 it is more important, 
however, to check if the change of the "narrowness", i.e. the ratio max l6 m i n — 1, to smaller 
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BHLUMI [nb] 


OBI+LMG [nb] 


SABSPV [nb] 


BHAGEN95 [nb] 




CAL02 Symmetric Narrow-Narrow 


.100 


95.458 ± .005 


95.259 ± .014 


95.363 ± .013 


95.287 ± .009 


.300 


95.233 ± .005 


95.048 ± .014 


95.157 ±.016 


95.065 ± .009 


.500 


94.841 ± .005 


94.672 ± .014 


94.792 ± .016 


94.680 ± .009 


.700 


93.520 ± .005 


93.347 ± .014 


93.513 ± .019 


93.354 ± .009 


.900 


87.359 ± .005 


86.899 ± .013 


87.396 ± .012 


86.958 ± .009 




CAL02 Asymmetric Narrow-Wide 


.100 


98.834 ± .003 


98.809 ± .010 


98.859 ± .017 


98.804 ± .009 


.300 


98.539 ± .003 


98.535 ± .010 


98.577 ±.017 


98.515 ± .009 


.500 


98.020 ± .003 


98.038 ± .010 


98.073 ± .019 


98.006 ± .009 


.700 


96.054 ± .003 


96.061 ±.010 


96.131 ± .018 


96.033 ± .009 


.900 


88.554 ± .003 


88.220 ± .009 


88.648 ± .015 


88.263 ± .009 
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Table 16: In this table/figure we show cross sections for various symmetric/asymmetric versions of the 
CAL02 ES, for matrix elements beyond first order. Z exchange, up-down interference and vacuum polar- 
ization are switched off. The center of mass energy is ^/s = 92.3 GeV. Not available x-sections are set to 
zero. The wide range is defined by 6\ w = 6\f ± 8 segm and 6i w = 2 f — S segm , and the narrow range by 



'1/ 



28 



segm 



and 6 2n = 02/ - 4£ segm ; 6, 



segm 



(02/ - 0i/)/16, 0i/ = 0.024 and 6 2f = 0.058 rad, 



respectively. 



values does not spoil the agreement of the table/figure 16. As we have already indicated, 
at LEP2 the decrease of the narrowness 8 max / ' 6 m in — 1 may cause a significant increase in 
the photonic radiative corrections. The relevant cross-check is done in table/figure 17. It 
represents the worst possible scenario at LEP2. The results are shown for the narrower version 
of the CAL02 ES, which we call CAL03, in the symmetric and asymmetric versions. As we 
see, BHLUMI and SABSPV differ again for the above ES by less than 0.2%. This result will 
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BHLUMI [nb] 


OBI+LMG [nb] 


SABSPV [nb] 


BHAGEN95 [nb] 




CAL02 Symmetric Narrow-Narrow LEP2 


.100 


8.088 ± .001 


8.052 ± .003 


8.074 ± .001 


8.058 ±.001 


.300 


8.074 ± .001 


8.039 ± .003 


8.061 ± .001 


8.044 ± .001 


.500 


8.048 ± .001 


8.014 ± .003 


8.039 ± .001 


8.018 ±.001 


.700 


7.989 ± .001 


7.954 ± .003 


7.986 ± .001 


7.958 ± .001 


.900 


7.574 ± .001 


7.515 ±.003 


7.582 ± .001 


7.522 ± .001 




CAL02 Asymmetric Narrow-Wide LEP2 


.100 


8.523 ± .001 


8.514 ± .002 


8.518 ±.001 


8.515 ±.001 


.300 


8.501 ± .001 


8.494 ± .002 


8.499 ± .001 


8.494 ± .001 


.500 


8.464 ± .001 


8.457 ± .002 


8.465 ± .001 


8.457 ± .001 


.700 


8.374 ± .001 


8.366 ± .002 


8.380 ± .001 


8.365 ±.001 


.900 


7.755 ± .001 


7.716 ±.002 


7.769 ± .001 


7.720 ± .001 
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Table 17: In this table/figure we show cross sections for for the symmetric/asymmetric CAL03 ES's (the 
narrower version of CAL02) for matrix elements beyond first order. Z exchange, up-down interference 
and vacuum polarization are switched off. The center of mass energy is ^/s = 176 GeV. Not available 
x-sections are set to zero. The wide range is defined by 6\ w = 6\f ± 6S segm and 6i w = 6\f ± 16S segm , and 



the narrow range by ln = 6 lf + 8S segm and 6 2n = #1/ ± I56 segm ; S segm 
and #2/ = 0.058 rad, respectively. 



(02/ - 0i/)/16, ^l/ = 0.024 



be used for estimating theoretical uncertainty of the SABH process at LEP2. Hybrid Monte 
Carlo's (OLDBIS ± LUMLOG^o) and BHAGEN95 are off of about 0.4% in the NN case but, 
noticeably, they are on the same ground as BHLUMI and SABSPV for the most interesting 
NW case. 
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BHLUMI [nb] 


SABSPV [nb] 


BHAGEN95 


VP+Z Bhlumi 




CAL02 Symmetric Wide-Wide 


.100 


136.975 ± .010 


136.831 ± .018 


136.861 ± .008 


5.140 ± .008 


.300 


136.576 ± .010 


136.453 ± .018 


136.482 ± .008 


5.126 ± .008 


500 


135.827 ± .010 


135.742 ± .018 


135.770 ± .008 


5.100 ± .008 


700 


132.962 ± .010 


132.928 ± .017 


132.927 ± .008 


4.994 ± .008 


Q00 


123.420 ± .009 


123.430 ± .018 


123.114 ± .009 


4.627 ± .008 




CAL02 Symmetric Narrow-Narrow 


1 00 


99.208 ± .009 


99.074 ± .017 


99.089 ± .011 


3.751 ± .007 


.300 


98.975 ± .009 


98.851 ± .021 


98.866 ± .011 


3.742 ± .007 


500 


98.570 ± .009 


98.477 ± .017 


98.479 ± .011 


3.728 ± .007 


700 

■ 1 uu 


97.198 ± .008 


97.147 ± .017 


97.128 ± .011 


3.678 ± .007 


Q00 


90.789 ± .008 


90.791 ± .016 


90.537 ± .011 


3.430 ± .007 




CAL02 Asymmetric Narrow- Wid 


e 


.100 


102.717 ± .006 


102.703 ± .017 


102.724 ± .010 


3.883 ± .004 


.300 


102.412 ± .006 


102.411 ± .017 


102.434 ± .010 


3.873 ± .004 


.500 


101.874 ± .006 


101.894 ± .017 


101.922 ± .010 


3.854 ± .004 


.700 


99.833 ± .005 


99.878 ± .017 


99.902 ± .010 


3.779 ± .004 


.900 


92.033 ± .005 


92.088 ± .016 


91.887 ± .011 


3.478 ± .004 



Table 18: Monte Carlo results for various symmetric/asymmetric versions of the CAL02 ES, for matrix 
elements beyond first order. Z exchange, up-down interference and vacuum polarization are switched ON. 
The center of mass energy is ^ = 92.3 GeV. Not available x-sections are set to zero. The wide range is 
defined by lw = lf + S segm and 2w = # 2 / - S segm , and the narrow range by ln = 6 lf + 2S segm and 
#2n = 02/ - ^segm) & segm = (#2/ - 0i/)/16, 0\j = 0.024 and 2 f = 0.058 rad, respectively. 

2.7.5 Z and vacuum polarization included 

In all the previous comparisons, the small contributions from s-channel Z-exchange and s- 
channel photon exchange diagrams were switched off in order to enhance the possibility of 
seeing more clearly the most important pure photonic higher order corrections. In the following 
part of numerical comparisons, we restore in the calculations the contributions from these 
5-channel Z-exchange and s-channel photon exchange diagrams, together with the effect of 
vacuum polarization. The comparison of various calculations is done for the semi-realistic ES 
CAL02 in the versions Wide- Wide, Narrow-Narrow and Narrow- Wide, as defined in Fig. 13. 
The resulting cross sections are shown for a LEP1 energy in Tab. 18 and Fig. 17. Again, 
BHLUMI and SABSPV, for values of the energy-cut variable in the experimentally interesting 
range 0.25 < z min < 0.75, agree within 0.1% for all the three versions of the ES (WW, NN and 
WN). BHAGEN95 is also in agreement, in this case, for all the three versions of the ES, due to 
a slightly bigger correction in these added contributions. We do not expect that switching on 
the small s-channel Z-exchange and s-channel photon exchange corrections would change our 
conclusions for LEP2. Vacuum polarization enters essentially only in the normalization of the 
SABH cross section, and Z contribution at LEP2 can be safely neglected. We therefore extend 
the validity of the above exercise to LEP2. 
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Figure 17: Monte Carlo results for various symmetric/asymmetric versions of the CAL02 ES, for matrix 
elements beyond first order. Z exchange, up-down interference and vacuum polarization are switched ON. 
The center of mass energy is ^ = 92.3 GeV. Not available x-sections are set to zero. In the plot, the 
0( a2 )exp S cross section c BHL from BHLUMI 4.x is used as a reference cross section. 

2.8 The total theoretical error for small-angle Bhabha scattering 



In this section we present some supplementary numerical material concerning higher order 
corrections from MC and non-MC programs, and we summarize on the total theoretical error 
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BHLUMI(alf2e) 


BHLUMI(alf3e) 


BHLUMI(alf2) 


NLLBHA(alf2) 


NLLBHA(alf3) 


NLLBHA(alf3p) 




(a) BARE1 


.100 


166.892 ± .006 


-.017 ±.000 


166.988 ± .021 


167.016 ±.017 


166.948 ± .000 


166.966 ± .000 


.300 


165.374 ± .006 


-.010 ±.000 


165.471 ± .021 


165.503 ± .017 


165.448 ± .000 


165.421 ± .000 


.500 


162.530 ± .006 


-.006 ±.000 


162.594 ± .021 


162.630 ± .016 


162.581 ± .000 


162.527 ± .000 


.700 


155.668 ± .006 


-.002 ±.000 


155.620 ± .020 


155.649 ± .015 


155.617 ±.000 


155.528 ± .000 


.900 


137.342 ± .006 


.004 ± .000 


137.191 ±.020 


137.205 ± .014 


137.201 ± .000 


137.063 ± .000 




(b) SICAL2 


.000 


132.816 ± .006 


-.017 ±.000 


132.912 ± .019 


.000 ± .000 


.000 ± .000 


.000 ± .000 


.200 


132.553 ± .006 


-.018 ±.000 


132.645 ± .019 


.000 ± .000 


.000 ± .000 


.000 ± .000 


.400 


131.985 ± .006 


-.019 ±.000 


132.061 ± .019 


.000 ± .000 


.000 ± .000 


.000 ± .000 


.600 


128.672 ± .006 


-.017 ±.000 


128.711 ±.019 


.000 ± .000 


.000 ± .000 


.000 ± .000 


.800 


119.013 ±.006 


-.012 ±.000 


119.014 ±.018 


.000 ± .000 


.000 ± .000 


.000 ± .000 
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Table 19: In this table/figure we show cross sections for LEP1 center of mass energy, y/s = 92 GeV. 
Results from BHLUMI and N LLBHA for the symmetric Wide-Wide ES's BARE1 and SICAL2 are shown. Not 
available x-sections are set to zero. In the table, column BHLUMI(alf2e) represents 0(a 2 ) exp BHLUMI 
4.02.a, col. BHLUMI(alf2) shows C(a 2 ) BHLUMI without exponentiation, col. BHLUMI(alf3e) shows 
missing 0(cl 3 )ll in BHLUMI 4. 02. a as calulated with the new (unpublished) version of LUMLOG, col. 
NLLBHA(alf2) shows 0(a 2 ) result from NLLBHA including NLL corrections, col. NLLBHA(alf3) is the 
previous plus 0(o. 3 )ll and col. N LLBHA(alf3p) is the previous plus light pair corrections. In the plot, the 
0(a 2 ) exp cross section cr KBF from BHLUMI 4. 02. a is used as a reference cross section (except for missing 
0(o. 3 )ll, for which we show (t/ctref)- 

for the SABH process at LEP1 and LEP2. 

Let us discuss again the size of the 0(a 3 L 3 ) and 0(a 2 L) corrections. In the next ta- 
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LEP1 


LEP2 


Type of correction/error 


Ref. [6] 


Present 


Present 


(a) Missing photonic 0(a 2 L) 


0.15% 


0.10% 


0.20% 


(a) Missing photonic 0(a 3 L 3 ) 


0.008% 


0.015% 


0.03% 


(c) Vacuum polarization 


0.05% 


0.04% 


0.10% 


(d) Light pairs 


0.04% 


0.03% 


0.05% 


(e) Z-exchange 


0.03% 


0.015% 


0.0% 


Total 


0.16% 


0.11% 


0.25% 



Table 20: Summary of the total (physical+technical) theoretical uncertainty for a typical calorimetric 
detector. For LEP1, the above estimate is valid for the angular range within 1° — 3°, and for LEP2 it covers 
energies up to 176 GeV, and angular range within 1° — 3° and 3° — 6° (see the text for further comments). 

ble/figure 19, we address this question showing once again some results from Tab. 14/Fig. 16, 
and adding some new numerical results from the BHLUMI event generator and the semianalyt- 
ical program NLLBHA for the unrealistic ES BARE1 and the realistic ES SICAL2, symmetric 
WW variants. First, let us recall that in Tab. 14/Fig. 16 the 0(a 3 L 3 ) effects were included 
through exponentiation in all calculations, but in most cases they were incomplete. In the case 
of BHLUMI, the recent version of LUMLOG 6 is able to answer the question: how big is the 
missing 0(a 3 L 3 ) in BHLUMI 4.02a. In table/figure 19 we see (black dots) that it is below 
0.01% for both BARE1 and SICAL2 ES's. According to our "scaling rules", we conclude that 
it is below 0.02% at LEP2. Hence, from the practical point of view, 0(a 3 L 3 ) in BHLUMI 
4.02a is complete. In table/figure 19 we also include, for the unrealistic BARE1 ES, numerical 
results from NLLBHA (stars), which includes complete 0(a 2 L) and 0(a 3 L 3 )ll corrections. 
The difference between BHLUMI (crosses) and NLLBHA (stars) should be, in principle, due to 
0(a 2 L) (and technical precision), because 0(a 3 L 3 ) should cancel completely. As we see, the 
above difference is within the "one per mil box", but for stronger cuts, z m i n = 0.9, it grows 
slightly beyond 0.1%. Luckily enough, we may push the above exercise in the interesting direc- 
tion - we have also in table/figure 19 the results from BHLUMI (circles) and NLLBHA (boxes), 
in which exponentiation and 0(a 3 L 3 )ll was removed completely. As we see, these results agree 
better, even for strong energy cut (z m i n = 0.9). Actually, this result (difference between boxes 
and circles) represents an interesting quantity: missing 0(a 2 L) in BHLUMI. The above result 
suggests that it is rather small, below 0.03%. One has to keep in mind that, if the above is 
true, then the former difference, with 0(a 3 L 3 )ll (crosses and stars), is a puzzle and needs to 
be examined further. In any case, the fact that all the four above results from BHLUMI and 
NLLBHA are within the "one per mil box" is interesting, encouraging and reinforcing our final 
conclusion that photonic corrections are under control within 0.1%. For the present time the 
above interesting comparison is limited to BARE1 ES. For SICAL2 and BARE1 ES's, we see 
that the difference between BHLUMI with and without exponentiation is quite sizeable, 0.08%, 
and from that we conclude that the inclusive Yennie-Frautschi-Suura exponentiation in BH- 

6 The new LUMLOG includes final state radiation (in addition to the initial) up to 0(a z L 3 )ll. It was 
discussed in the Bhabha Working Group and will be included in the next release of BHLUMI. 
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LUMI is necessary and instrumental for getting good control over the 0(a 3 L 3 )ll corrections, 
even if they are not complete in the matrix element. As a matter of fact, all the other MC codes 
involved in the present study include exponentiation, and so are on a firm ground from this 
point of view. In table/figure 19, we show also results from NLLBHA including pair production 
in addition to the 0(a 2 L) and 0(a 3 L 3 ) corrections ("plus" marks in the plot). The difference 
between pluses and stars represents the net effect of the light fermion pair production. For the 
BARE1 ES, with z m i n in the experimentally interesting range, it is 0.06% at most. We expect 
this effect to be about a factor of two smaller for calorimetric ES's. 

The total theoretical error for the SABH process at LEP1/LEP2 is summarized in table 20. 
The errors in the table are understood to be with respect to the cross section calculated for any 
typical (asymmetric) ES, for the LEP1 experiment in the angular range 1° — 3°, with respect to 
the cross section calculated using BHLUMI 4.02a. In the case of LEP2, the estimate extends 
to the angular range 3° — 6°, and to the case of the angular range about twice narrower than 
usual (see the discussion of the numerical results in the previous sections). The entries include 
combined technical and physical precision. In this table, entry (a) for Missing 0(a 2 L) is based 
mainly on the agreement between BHLUMI and SABSPV, as seen in tables 14, 16 and 18. 
It should be stressed that we rely on the agreement between BHLUMI and SABSPV for all 
the three types of ES, Wide- Wide, Wide-Narrow and Narrow Narrow. The agreement between 
BHLUMI and SABSPV is now better than the one between BHLUMI and OLDBIS+LUMLOG 
used in the previous best error estimate of Ref . [6] . Noticeably, albeit the agreement between 
BHLUMI on the one side and BHAGEN95/(OLDBIS + LUMLOG) on the other side is not 
always below 0.1% for all the ES's considered, it is at least for the experimentally most interest- 
ing NW case. This fact is a further reinforcement of the present theoretical error estimate for 
the SABH process in the NW case, and it is a suggestion for the experimentalists to continue 
to choose the NW-ES's. The fact that for the unrealistic ES BARE1 the difference between 
BHLUMI and NLLBHA, see fig. 19, is also within 0.1% confirms this evaluation. Entry (b) is 
based on table/figure 19. In entry (c), the new improved uncertainty of the vacuum polariza- 
tion is taken from Tab. 7. We take the biggest of the results from refs. [33,34]. The light pair 
production uncertainty, entry (c), is based on new estimates reported during the workshop (see 
Ref. [7,8,12,15] and Ref. [26,28]; see also table/figure 19). In tab. 20, we quote for LEP1 the 
present error due to light fermion pairs contribution to be 0.03%. This is based on all the refer- 
ences quoted above and on the discussion during the WG meetings [29]. The previous estimate 
in Ref. [6] is therefore confirmed and improved slightly. This is under the assumption that the 
pair effect is corrected for at least in the LL approximation. If the effect is not corrected for 7 , 
then we recommend to use for LEP1 0.04% as an estimate for the missing pair effect (0.06% 
for LEP2). The material presented at the workshop suggests that the final uncertainty of the 
light pair contribution will be at the level of 0.015%. In entry (e), the reduced uncertainty of 
the Z-exchange contribution is based on Ref. [59], work done during this Workshop. 

The improvement of the theoretical luminosity error from 0.16% down to 0.11% is basically 

7 Production of the light pairs is not included in the standard version of BHLUMI. It is implemented only in 
the testing unpublished version [29]. 
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due to successful comparisons of the programs BHLUMI and SABSPV for a wide range (WW, 
NN and NW) of experimentally realistic ES's (SICAL2), and also due to an encouraging (al- 
though limited to the unrealistic ES BARE1) comparison of un-exponentiated BHLUMI and 
NLLBHA in table/figure 19. Furthermore, the agreement of BHLUMI, SABSPV, BHAGEN95 
and (OLDBIS+LUMLOG) within that same 0.11% error in the NW-ES recommends safely this 
choice in the experimentally relevant cases. At last, the analysis described in subsection 2.6 
shows that the actual Bhabha selections used by the LEP experiments to measure the acceler- 
ator luminosity minimize the sensitivity to 0(a 2 ) radiative corrections, thus putting the above 
conclusions on an even firmer ground. We would like to stress very strongly that the above 
new estimate 0.11% of the total luminosity error is based on new results which, although pretty 
stable numerically, are generally still quite fresh and they are unpublished. We expect these new 
results to be published in journals shortly after the workshop, together with the corresponding 
computer programs. 

The total theoretical error for the SABH process at LEP2 is also summarized in Tab. 20. 
We assume that the cross section is calculated for any typical (asymmetric) ES for LEP2 
experiment, in the typical angular range 1° — 3° or 3° — 6°. The error estimate covers also 
the "worst case scenario" of the super-narrow angular range (see the example of ES CAL03 
in table/figure 17). In entry (a), the estimate of the total photonic uncertainty is based again 
upon the agreement between BHLUMI and SABSPV on all the variants of ES's considered, 
and reinforced by the fact that BHAGEN95/OLDBIS+LUMLOG are on the same ground as 
BHLUMI and SABSPV in the experimentally more interesting NW case (see tables /figures 15 
and 17). Note that, sometimes, in the case of other angular range 3° — 6° and higher energies, 
the "scaling laws" from the introduction were used instead of direct calculation to extend 
the actual numerical results to these situations (see the comments accompanying the relevant 
tables/plots). We do not see much danger in this because, usually, the large safety margin close 
to a factor of two was present. Entry (b) is produced out of LEP1 result using the "scaling 
rule". The vacuum polarization for LEP2 case in the Tab. 20 is taken from Tab. 5 at the 
\t\ = 36 GeV 2 , corresponding to LEP2 energy and the angle of 6 m in = 60 mrad. 



Type of correction/error 


Error estimate 


(a) Missing 0(a 2 L), 0(a 3 L 3 ) 


< 0.010 % 


(b) Technical precision (photonic) 


0.040% 


(c) Vacuum polarization 


0.030% 


(d) Light fermion pairs 


0.015% 


(e) Z-contribution 


0.010% 


Total 


0.053% 



Table 21: Future projection of the total (physical+technical) theoretical uncertainty for a typical calori- 
metric detector, within the 1° — 3° angular range at LEP1 energies. 

Finally, in view of all the work reviewed during the workshop, we are also able to estimate 
the precision which will be attained in the next step. It is shown schematically in table 21. 
At the time when Monte Carlo programs will include the matrix element from 0(a 2 L), the 
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uncertainty due to higher order corrections will be negligible. The dominant contribution will 
be of technical origin and we think that, as we have seen from 0(a 1 ) comparisons, it can 
be reduced to 0.04% (provided we can successfully tune two independent Monte Carlo event 
generators at that precision level, for the same or very similar 0(a 2 ) matrix elements). The 
vacuum polarization is now taken according to Ref. [33], and from the discussions during the 
workshop meetings it was obvious that a further reduction of the uncertainty due to pairs and 
Z-exchange is also possible. The corresponding work is in progress. 



3 Large-angle Bhabha scattering 

In the present section the LABH process is considered, both at LEP1 and LEP2. The aim of 
the study, rather than updating the conclusions of Ref. [1] concerning the theoretical accuracy 
of the LABH process at LEP1, is twofold: on the one hand, the comparison between the semi- 
analytical benchmarks and the Monte Carlo codes used by the LEP collaborations; on the other 
one, the study of the LABH process at LEP2, accompanied by the development of dedicated 
software. 



3.1 Physics 

The main physics interest of Bhabha scattering measurements at large angles (say 8 > 40°) 
around the Z resonance is a precise test of the electroweak sector of the Standard Model. In this 
angular region more than 80 % of the cross section is due to resonant s-channel Z exchange. For 
y/s = Mz the interference contributions between s-channel Z exchange and the other diagrams 
either vanish or are completely irrelevant, and the s-channel photon exchange contribution is 
small (~ 5 X 10~ 3 of the Z exchange cross section). The only other relevant contribution is 
^-channel photon exchange. For electroweak analyzes, one thus subtracts the ^-channel and 
s — t interference contributions from the large-angle experimental data, typically calculated 
using the ALIBABA [60] semi- analytical (SA) program. After correcting for the effects of real 
and virtual photon radiation using the analytical programs MIBA [61,62], TOPAZ0 [63,64] or 
ZFITTER [65,66], the Z exchange cross section a z may be extracted. For = Mz, cr% = <^ ak 
where: 

■peak _ 127Tr e 



Mfr| 



(2) 



For the other charged lepton pair decay modes, fi + fi~, t + t~ , of the Z the quantity Y 2 e in Eqn. (2) 
is replaced by T e T^, r e r T , respectively, while for hadronic (qq) decays it is replaced by T e Th a d- 
Thus the electronic width of the Z, T e , which appears in the cross section for all decay modes of 
the Z, is measured directly and with improved sensitivity (because in this case a^ ah oc Y\~) only 
in large-angle Bhabha scattering. It is worth noting, however, that in principle the so called 
^-channel subtraction is not unavoidable. Actually, the program TOPAZ0 [63, 64] could be 
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used to fit directly the data for large-angle Bhabha scattering without relying upon ^-channel 
subtracted data. 



The resulting sensitivity of the backward-forward charge asymmetry in large-angle Bhabha 
scattering to the important electroweak parameters 8 AK top and Ak higgs 

A top _ 3a/2G^c| 1 2 

AK ~ 16*' ^ « {6} 

higgs V2G»M* Wf 10, fM H \ 2 5, 



An = ^L(_ f --) (4) 

is similar to that of the other dilepton channels fi + fi~ , t + t~. The above formulae of course 
indicate only the leading dependence of the one-loop corrections on the masses of the top-quark 
and the Higgs boson. Actually, at the nowadays precision level a complete electroweak library 
is mandatory [1]. 

At the Z peak the purely QED corrections to the large-angle Bhabha cross section are, for 
typical experimental cuts [67]: 0(a), -30 % ; (9(a 2 ), +4 %. These corrections are much larger 
than those in small-angle Bhabha scattering when typical 'wide'/'narrow' cuts are used [32]: 
(9(a), +5 % ; (9(a 2 ), -1.4 %. Thus theoretical errors on QED radiatively corrected cross sections 
are expected to be considerably larger in large-angle than in small-angle Bhabha scattering. 
This is indeed found to be the case in the comparisons between different codes shown below. 

In the energy regime of LEP2, the Z-boson effects on the large-angle Bhabha cross section 
are much smaller than at LEP1. Actually, before entering the details of the comparisons, 
it is worth noting that large-angle Bhabha scattering shows very different physical features 
depending on the energy regime at which it is considered. As can be seen from Fig. 18, around 
the Z peak the cross section is largely dominated by Z-boson annihilation, whereas, already 
some GeV off resonance, the cross section is largely dominated by ^-channel photon exchange. 
From this point of view, large-angle Bhabha scattering at LEP2 is much more similar to small- 
angle Bhabha scattering than to large-angle Bhabha scattering at LEP1. Hence, at LEP2 the 
large-angle Bhabha cross section cannot be a useful tool for precise tests of the electroweak 
sector of the Standard Model, but rather for general QED tests. 

The state-of-the-art of large- angle Bhabha scattering up to now can be found in Ref. [1]. In 
that paper an extensive comparison between two semi- analytical codes, namely ALIBABA [60] 
and TOPAZ0 [63,64], is shown. On the other hand, although extensive in the sense that cross 
sections and asymmetries are considered, that comparison is in some sense limited: actually it 
involves only semi- analytical codes, on very simple, academic ES's, only at the Z peak. 

In view of the above considerations, the tasks of the present Working Group, as far as 
large-angle Bhabha scattering is concerned, are the following ones: 

• involving in the comparisons also the Monte Carlo (MC) codes today available and used 
by the LEP collaborations; 



8 See [1] for a discussion of pseudo-observables for precision calculations at the Z peak. 
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Figure 18: The relative contributions to the integrated cross section at the Born level. The individual 
contributions are, from top to bottom on the right-hand side of the plot: f(t)f(t), Z(t)Z(t), 7(5)7(5), 
Z{s)Z{s), Z{s)i{s), Z{s)Z{t), 1 {s)Z{t), Z{s)i{t), f(t)Z(t), and «f(s)<y(t). 

• considering also more realistic, albeit simple, ES's; 

• providing results also for the LEP2 energy range, eventually developing dedicated soft- 
ware. 

The ES's considered in the present study are the following ones: 



• BARE - This ES, for the sake of simplicity, is defined exactly as in [1], namely 40° < 
< 140°, 0° < tf+ < 180°, $™ c a oU = 10°, 25° and E min = 1 GeV for both electron and 
positron; 



• CALO - This ES is defined as above, but with E m i n = 20 GeV for the final fermion energy, 
which is the electron(positron) energy if there are no photons nearby, whereas it is the 
electron(positron) plus photon energy if the photon is within a cone of semi-aperture 1° 
from the electron(positron). 



For all the cases considered, the input parameters are Mz = 91.1887 GeV, m t = 174 GeV, 
m H = 300 GeV and a s {M z ) = 0.124. The predictions by ALIBABA are taken from Ref. [1]. 
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Let us now briefly summarize the features of the codes involved in the study. Here only 
the general features will be highlighted; for more details the reader is referred to the original 
literature or to the write-ups presented at the end of this section. 

ALIBABA [60] - It is a semi- analytical code, implementing exact 0(a) QED and weak cor- 
rections. The higher-order QED corrections consist of leading log 0(ct 2 ) corrections plus soft- 
photon exponentiation. Moreover, the weak 0(a) corrections are folded with the leading log 
structure functions. The matching between the exact 0(a) QED matrix element and the 
higher order corrections is performed in additive form. The electroweak library is not up to 
date. Nonetheless, the code has to be considered as a benchmark. 

BHAGEN95 [43] - It is a Monte Carlo integrator for both small and large-angle Bhabha scat- 
tering. The value for the cross section is obtained from the event generator BHAGEN94, 
a structure function based program for all orders resummation, including complete photonic 
0(a) and leading logarithmic 0(a 2 L 2 ) corrections in all channels, and all relevant electroweak 
corrections according to BHM/WOH basic formulae from Ref. [1]. The approximations, intro- 
duced with the collinear kinematics of initial and final radiation and in its angular distribution, 
are eliminated for the one hard photon emission by substitution with the exact calculation. 

BHAGENE3 [67, 68] - It is a Monte Carlo event generator for large-angle Bhabha scattering 
and muon pair production. The program includes one-loop and the most important two- 
loop electroweak as well as QED radiative corrections. The 0(a) QED correction uses the 
exact matrix element. Higher order QED corrections are included in an improved soft photon 
approximation with exponentiation of initial state radiation. Up to three hard final state 
photons are generated. Events are generated in the full final state phase space including explicit 
mass effects in the region of collinear mass singularities. The minimum scattering angle for 
percent level cross section accuracy is 10°. Extensive use is made in the program of one and 
two dimensional look-up tables for fast, flexible and efficient Monte Carlo generation. The 
program was designed for the Z peak region but may also be used at LEP2 energies. 

BHWIDE [69] - It is a new Monte Carlo event generator for large-angle Bhabha scattering 
at LEP1/SLC and LEP2. It includes multiphoton radiation in the framework of 0(a) YFS 
exponentiation. The 0(a) virtual (both weak and QED) corrections are in the current version 
taken from ALIBABA. The program provides the full event in terms of particle flavors and 
their four-momenta with an arbitrary number of radiative photons. In many aspects it is 
similar to the program BHLUMI for small-angle Bhabha scattering and can be considered as 
its extension to large angles. It has been checked that for the pure QED process BHWIDE at 
O(a) (no exponentiation) agrees with the MC program OLDBIS within a statistical accuracy 
of 0.05%. 

SABSPV [46] - It is a new Monte Carlo integrator, originally designed for small-angle Bhabha 
scattering, but adapted to the treatment of large-angle Bhabha scattering at the LEP2 energy 
range. It is based on a proper matching of the 0(a) corrected cross section for ^-channel 
photon exchange and of the leading logarithmic results in the structure function approach. 
The matching is performed in a factorized form, in order to preserve the classical limit. At 
present, the effect of up-down interference in the *y(t) — *y(t) contribution is not taken into 
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account and all the other contributions are corrected at the leading logarithmic level. Due to 
the present approximations, the theoretical accuracy of the code is of the order of 1%, as far 
as large-angle Bhabha scattering at LEP2 is concerned. 

TOPAZO [63,64] - It is a semi- analytical code, developed for precision physics at LEP1. It 
includes the state-of-the-art concerning weak and QCD corrections, according to Ref. [1]. As 
far as QED corrections are concerned, they are exactly treated at 0(a) for s-channel processes 
(leptonic and hadronic), at the leading logarithmic level for pure ^-channel and s-t contributions 
in the Bhabha scattering case. On top of this, higher order QED corrections are taken into 
account in the structure functions approach, in a factorized form in order to preserve the clas- 
sical limit. A particular effort has been performed in order to implement as much analytically 
as possible the experimental cuts typically applied by the LEP collaborations. 

UNIBAB [70] - It is a full Monte Carlo event generator that was originally designed for large- 
angle Bhabha scattering at LEP1 and SLC energies. The QED corrections are implemented in a 
fully factorized form by assuming s-channel dominance and using photon shower algorithms for 
initial- and final-state radiation, and therefore exponentiation of soft photons and resummation 
of the logarithms from multiple emission of hard collinear photons is automatic. QED initial- 
final interference corrections are not yet implemented. The electroweak corrections are based 
on a library also used by ALIBABA, but updated to include the leading ra^-dependence and 
higher order QCD corrections to the Z width. 

3.2 On Z peak (LEP1) 

The situation of the comparisons for LEP1 is summarized in Figs. 19 (BARE) and 20 (CALO) 
and corresponding tables. Conventionally, the reference cross section with respect to which the 
relative deviations are computed is taken from TOPAZO. It has to be stressed that this choice 
has no particular meaning at all. 

Let us begin with commenting the situation of Fig. 19, i.e. for the BARE ES. As far as 
the comparison between the two semi- analytical codes, ALIBABA and TOPAZO, is concerned, 
the agreement is better than 0.1% at the Z peak (energy points n. 4 and 5, corresponding to 
the smallest experimental error, which is of the order of 0.3% statistical and 0.3% systematic), 
and deteriorates on the wings, where, on the other hand, the experimental error is larger (for 
instance, at peak±2 GeV the experimental error is of the order of 1% statistical and 0.3% 
systematic). Note that the worst situation is for maximum acollinearity cut of 10°, above the 
Z peak, where the codes differ from one another of about 1%: this difference is due to higher 
order QED effects, as pointed out in Ref. [71] (factorized versus additive formulation). As far 
as the Monte Carlo codes BHAGENE3 and BHWIDE are concerned, their agreement with the 
semi- analytical codes at peak is within few per mil, whereas off peak BHWIDE is within 1% 
and BHAGENE3 can deviate up to 2%. 

The situation for the more realistic case, the CALO ES (Fig. 20), is generally better from 
the point of view of the SA/MC comparisons. Note that ALIBABA is no more involved, since 
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TOPAZ0 
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BHAGEN95 




(a) BARE acoZ raa3! = 10° 


1. 


88.45 


.4579 ± .0003 


.4560 ± .0004 


.4495 ± .0016 


.4575 ± .0003 


.4556 ± .0002 


A. 


89.45 


.6452 ± .0002 


.6429 ± .0006 


.6334 ± .0023 


.6440 ± .0003 


.6403 ± .0003 


Q 
0. 


90.20 


.9115 ± .0002 


.9087 ± .0008 


.8997 ± .0033 


.9090 ± .0004 


.9026 ± .0004 


4 


91.19 


1.1846 ± .0002 


1.1797 ± .0010 


1.1847 ± .0033 


1.1840 ± .0004 


1.1715 ± .0005 


c. 
o . 


91.30 


1.1639 ± .0002 


1.1592 ± .0009 


1.1667 ± .0033 


1.1636 ± .0005 


1.1514 ± -0005 


a 
u. 


91.95 


.8738 ± -0002 


.8711 ± .0007 


.8856 ± .0028 


.8769 ± .0003 


.8664 ± .0003 


7 


93.00 


.4771 ± -0002 


.4761 ± .0005 


.4808 ± .0019 


.4814 ± .0001 


.4756 ± .0002 


Q 
O. 


93.70 


.3521 ± -0002 


.3512 ± .0004 


.3521 ± .0013 


.3556 ± .0001 


.3522 ± .0001 




(b) BA] 


ElE acolrnaz, — 25° 


1 

1. 


88.45 


.4854 ± -0003 


.4808 ± .0005 


.4699 ± .0016 


.4833 ± .0003 


.4833 ± .0003 


o 
z. 


89.45 


.6746 ± -0003 


.6699 ± .0006 


.6593 ± .0023 


.6727 ± .0003 


.6727 ± .0003 


Q 
G. 


90.20 


.9438 ± -0003 


.9387 ± .0008 


.9279 ± .0033 


.9425 ± .0003 


.9425 ± .0003 


4. 


91.19 


1.2198 ± .0003 


1.2130 ± .0010 


1.2169 ± .0034 


1.2187 ± .0004 


1.2187 ± .0004 


5. 


91.30 


1.1989 ± .0003 


1.1924 ± .0010 


1.1995 ± .0034 


1.1982 ± .0004 


1.1982 ± .0004 


6. 


91.95 


.9054 ± .0002 


.9011 ± .0007 


.9124 ± .0026 


.9089 ± .0003 


.9089 ± .0003 


7. 


93.00 


.5040 ± .0002 


.5013 ± .0005 


.4996 ± .0019 


.5054 ± .0002 


.5054 ± .0002 


8. 


93.70 


.3777 ± .0002 


.3749 ± .0004 


.3689 ± .0013 


.3782 ± .0001 


.3782 ± .0001 



BARE acoL 



10° 



BARE acoL 



25° 



-.02 - 



<7 ~ |7 REF 
CT REF 



I ' ' 1 1 I ' ' 1 1 I ' 1 1 
O BHWIDE 

★ BHAGENE3 
o ALIBABA 
n BHAGEN95 
X TOPAZ0<r REF 




.02 



.01 



.01 



--.02 - 



g - g RF,F 
°REF 



I 1 ' 1 1 I 1 1 1 1 I 1 1 1 
O BHWIDE 

★ BHAGENE3 
o ALIBABA 
a BHAGEN95 
X TOPAZ0<7 REF 




5. 6. 7. 8. 



Figure 19: Monte Carlo results for the BARE ES, for two values (10° and 25°) of acollinearity cut. Center 
of mass energies (in GeV) close to Z peak. In the plots, the cross section <tref from TOPAZ0 is used as 
a reference cross section. Cross sections in nb. 



it cannot manage calorimetric measurements, whereas UNIBAB appears (it is slow for very 
small minimum fermion energy and therefore it did not contribute to the BARE case). On 
peak, the agreement between the codes is at the few per mil level; off peak BHWIDE is within 



48 



IN O. 


Ecm 


TOPAZ0 


BHWIDE 


BHAGENE3 


UNIBAB 


BHAGEN95 




(a) CALO acoZ raa3! = 10° 


1. 


88.45 


.4533 ± .0004 


.4523 ± .0004 


.4467 ± .0008 


.4490 ± .0010 


.4501 ± .0002 


A. 


89.45 


.6387 ± .0004 


.6377 ± .0006 


.6302 ± .0011 


.6358 ± .0012 


.6326 ± .0003 


Q 
0. 


90.20 


.9023 ± .0003 


.9016 ± .0008 


.8920 ± .0015 


.9021 ± .0014 


.8918 ± .0004 


4 


91.19 


1.1725 ± .0001 


1.1707 ± .0010 


1.1767 ± .0021 


1.1772 ± .0016 


1.1582 ± .0005 


c. 
o . 


91.30 


1.1520 ± .0001 


1.1505 ± .0009 


1.1571 ± .0020 


1.1559 ± .0016 


1.1385 ± -0005 


a 
u. 


91.95 


.8649 ± -0001 


.8646 ± .0007 


.8795 ± .0015 


.8689 ± .0012 


.8579 ± .0003 


7 


93.00 


.4723 ± .0001 


.4725 ± .0005 


.4796 ± .0008 


.4733 ± .0008 


.4719 ± .0002 


Q 
O. 


93.70 


.3486 ± -0001 


.3486 ± .0004 


.3507 ± .0006 


.3486 ± .0007 


.3498 ± .0001 




(b) CA1 


LO acolrnaz, — 25° 


1 

1. 


88.45 


.4769 ± -0004 


.4742 ± .0004 


.4696 ± .0008 


.4733 ± .0010 


.4717 ± .0002 


o 
z. 


89.45 


.6638 ± .0003 


.6615 ± .0006 


.6556 ± .0011 


.6619 ± .0012 


.6554 ± .0003 


Q 
G. 


90.20 


.9297 ± -0003 


.9278 ± .0008 


.9207 ± .0012 


.9302 ± .0014 


.9164 ± .0004 


4. 


91.19 


1.2025 ± .0003 


1.1994 ± .0010 


1.2074 ± .0021 


1.2073 ± .0016 


1.1845 ± .0005 


5. 


91.30 


1.1819 ± .0003 


1.1790 ± .0010 


1.1879 ± .0021 


1.1860 ± .0016 


1.1647 ± .0005 


6. 


91.95 


.8924 ± .0003 


.8909 ± .0007 


.9058 ± .0016 


.8965 ± .0012 


.8817 ± .0003 


7. 


93.00 


.4964 ± .0003 


.4954 ± .0005 


.5004 ± .0009 


.4976 ± .0008 


.4929 ± .0002 


8. 


93.70 


.3717 ± .0003 


.3704 ± .0004 


.3690 ± .0006 


.3720 ± .0007 


.3701 ± .0001 



-.02 - 



CALO acoL 



10° 



CALO acoL 



25° 



<7 ~ |7 REF 
CT REF 



I ' ' 1 1 I ' ' 1 1 I ' 1 1 
O BHWIDE 

★ BHAGENE3 
o UNIBAB 
n BHAGEN95 
X TOPAZ0<r REF 



+ 



+ 




.02 



.01 



.01 



--.02 - 



r ~ g RF,F 
CT REF 



I ' ' ' 1 I 1 ' 1 1 I 1 1 1 1 I 1 1 1 1 I 1 
O BHWIDE 

★ BHAGENE3 

o UNIBAB 

n BHAGEN95 ^ 



X TOPAZOcr 



REF 




Figure 20: Monte Carlo results for the CALO ES, for two values (10° and 25°) of acollinearity cut. Center 
of mass energies (in GeV) close to Z peak. In the plots, the cross section <tref from TOPAZ0 is used as 
a reference cross section. Cross sections in nb. 



0.5% from TOPAZ0, whereas UNIBAB deviates up to 1% below peak, and BHAGENE3 can 
differ from TOPAZ0 by about 2%. 
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BHAGEN95 is within 1.5% everywhere for both the BARE and CALO ES's around the 
Z peak. The agreement is better few GeV above and below the Z resonance. However the 
implementation of the complete weak and QCD library is very recent and still under tests. 



3.3 Far off Z peak (LEP2) 



No. 


BHWIDE 


TOPAZ0 BHAGENE3 


UNIBAB 


SABSPV 


BHAGEN95 




(a) CALO acol max = 10° 


1. 
2. 
3. 


35.257 ± .040 
29.899 ± .034 
25.593 ± .029 


35.455 ± .024 
30.024 ± .020 
25.738 ± .015 


34.690 ± .210 
28.780 ± .170 
24.690 ± .150 


34.498 ± .157 
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24.976 ± .115 


35.740 ± .080 
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25.960 ± .060 


35.847 ± .022 
30.352 ± .017 
26.007 ±.014 




(b) CALO acol max = 25° 


1. 
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3. 
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28.929 ± .036 
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27.840 ± .160 


39.521 ±.158 
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34.331 ± .020 
29.437 ± .015 
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Figure 21: Monte Carlo results for the CALO ES, for two values (10° and 25°) of acollinearity cut. Center 
of mass energies close to W^-pair production threshold (Ecm- 1- 175 GeV, 2. 190 GeV, 3. 205 GeV). In 
the plots, the cross section <tref from BHWIDE is used as a reference cross section. Cross sections in pb. 

The situation of the comparisons for LEP2 is shown in Fig. 21 (CALO) and corresponding table. 
Conventionally, the reference cross section with respect to which the relative deviations are 
computed is taken from BHWIDE. It has to be stressed again that this choice has no particular 
meaning at all. Note that TOPAZO has been developed in the Z-dominance approximation, and 
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UNIBAB does not include initial-final interference effects, so that their results are at the leading 
logarithmic level in the LEP2 energy range. A new entry appears, namely SABSPV, which has 
been conceived for small-angle Bhabha scattering, and further improved for large-angle Bhabha 
at LEP2. 

BHAGEN95, BHWIDE and SABSPV stay within 2% from one another. More precisely, 
SABSPV is steadily around 1.5% above BHWIDE and 0.5% below BHAGEN95. BHAGENE3, 
for both the acollinearity cuts considered, can deviate from the reference cross section up to 
5%. 

TOPAZ0 and UNIBAB show deviations from the reference cross section (up to 2% above and 
3% below, respectively) which depend on the acollinearity cut, and can be presumably traced 
back to the approximations intrinsic in these Z-peak designed codes. Anyway, the deviations of 
the two codes from the reference cross section are consistent with what can be expected from 
leading logarithmic results. 

4 Short-write-up's of the programs 

The aim of the following short- write-up's is to provide quick reference for the reader on basic 
properties of all event generators used in the numerical comparisons throughout this article. The 
intention was that details are given only on new and/or unpublished features of the programs 
(including bugs) while other features are described in general terms with help of references to 
published works. 

4.1 BHAGEN95 

AUTHORS: 

M. Caffo INFN and Dipartimento di Fisica dell'Universita, Bologna, Italy 

caf f o@bo . inf n . it 

H. Czyz University of Silesia, Katowice, Poland, INFN, Universita, Bologna, Italy 

czyzObo . inf n . it 

E. Remiddi INFN and Dipartimento di Fisica dell'Universita, Bologna, Italy 
remiddiObo . inf n . it 

GENERAL DESCRIPTION: 

BHAGEN95 is a collection of three programs to calculate the cross-section for Bhabha scattering 
for small and large scattering angles at LEP1 and LEP2 energies. In its present form the 
integrated cross-section a for a given selection of cuts is calculated as 

a = <r(BHAGEN94) - <r ff (BH94-FO) + ^(BHAGEN-IPH) . (5) 

<t(BHAGEN94) is the integrated cross-section obtained with the Monte Carlo event generator 
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BHAGEN94 [2,72-76], a structure function based program for all orders resummation, includ- 
ing complete photonic O(ct) and leading logarithmic 0(ct 2 L 2 ) corrections in all channels, and 
all relevant electroweak corrections according to BHM/WOH basic formulae from [1]. Approx- 
imations are introduced with the collinear kinematics of initial and final radiation and in its 
angular distribution. 

cr ff (BH94-FO) is the integrated cross-section of O(a) for one hard photon emission obtained 
with the Monte Carlo event generator BH94-FO [76], the O(a) expansion of BHAGEN94. 
cr^BHAGEN-lPH) is the integrated cross-section obtained with the one hard photon com- 
plete matrix element and exact kinematics, implemented in the Monte Carlo event generator 
BHAGEN-1PH [77]. 

The subtraction of 0- H (BH94-FO) and its substitution with ^(BHAGEN-IPH) is to eliminate 
the error in the contribution coming from the one hard photon emission. 

FEATURES OF THE PROGRAM: 

The three programs provide cross-sections, which are summed as in Eq. (5) or used to obtain 
other quantities, such as forward-backward asymmetry. Due to the mentioned substitution 
procedure, the event generator feature of the constituent programs can not be profited and the 
use is simply that of a Monte Carlo integrator. 

At small-angle we estimate the accuracy in the cross-section evaluation, which comes from the 
uncontrolled higher orders terms 0(a 2 L) and 0(a 3 L 3 ) and from the incertitude in 0(a 2 L 2 ) 
s — t interference to amount comprehensively to a 0.1%. The error due to approximate two 
hard photon contribution (strongly dependent on the imposed cuts) is estimated on the basis 
of the correction required for the one hard photon contribution times (3(s) = 0.1, to account 
for the increase in perturbative order. All included we estimate at small-angle an accuracy of 
the order of 0.15% for loose cuts (z m i n = 0.3) and of 0.45% for sharp cuts (z m i n = 0.9) for both 
LEP1 and LEP2 energies. 

At large-angle we estimate the 0(a 2 L 2 ) s — t interference accuracy up to 1% (depending on 
cuts) at LEP1 energy, but much smaller at LEP2. The error coming from the approximate 
treatment of two hard photon emission is estimated as above and is smaller for more stringent 
acollinearity cut. All included we estimate an accuracy of the order of 1% for both LEP1 and 
LEP2 energies. 

HOW DOES THE CODE WORK: 

The three programs run separately. They provide initialization and fiducial volume definition 
according to input parameters, then starts the generation of events according to some variables 
which smooth the cross-section behavior. Rejection is performed through the routine TRIGGER, 
where the special cuts can be implemented. The programs stop when the requested number of 
accepted events is reached or alternatively when the requested accuracy is obtained. 

INPUT CARD: 

The following data have to be provided in input: mass of the Z°, mass of the top quark, 
mass of the Higgs, value of as(M§), value of Tz, the b earn energy E\, eam ^ the minimum energy 
for leptons E m i n (larger than 1 GeV), minimum and maximum angle for the scattered electron 
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(positron) with the initial electron (positron) direction, maximum acollinearity allowed between 
final electron and positron, number of accepted events to be produced, numbers to initialize the 
random number generator. One may also switch on or off: pairs production, the channels to 
be considered and the recording of the events. For O(a) programs one has also to specify the 
minimum and maximum energy allowed for the photon. For the input of BHAGEN-1PH one 
has to give also the maximum acoplanarity, and minimum angles of the emitted photon with 
initial and final fermion directions, if one wants to exclude the contributions with the collinear 
photons. 

DESCRIPTION OF THE OUTPUT: 

Each program return the input parameters and the values of the cross-section obtained with 
weighted and unweighted events, with the relative statistical variance (one standard deviation). 
Of course due to the efficiency the weighted cross-section is usually much more precise than 
the unweighted one. The total integrated cross-section is then calculated according to Eq. (5). 

AVAILABILITY: 

On request to the authors and to be posted on WWW at http://www. bo.infn.it. / 

4.2 BHAGENE3 

AUTHORS: 

J. Field Departement de Physique Nucl. et Corpusculaire Univ. de Geneve 

jfield@cernvm.cern . ch 
T. Riemann DESY, Platanenallee 6, D-15738 Zeuthen 

riemannOif h . de 

GENERAL DESCRIPTION: 

BHAGENE3 is a Monte Carlo event generator for muon pairs (at all angles) or for Bhabha 
scattering in the large angle region (8 > 10°). The program, which is intended for use at, 
or above, the Z peak region contains all tree level diagrams with complete one loop and the 
leading two loop virtual corrections [78-81] The running a is included with the correct scale in all 
amplitudes. The 0(a) QED correction uses the exact ZZ7 matrix elements [82,83]. Higher order 
QED corrections are included in an improved soft photon approximation with exponentiation 
of initial state radiation. Events with up to three hard photons are generated in the full 
kinematically allowed phase space including explicit mass effects for near collinear photon 
radiation. If n^, are the respective numbers of initial and final state photons, the different 
final state topologies generated are: n^/n^ = 0/0, 1/0, 0/1, 2/0, 1/1, 0/2, 0/3 . Initial/Final 
state interference effects are taken into account only to 0(a). The photon energies are described 
by scaled variables: = 2El f /y /f s < 1. For < y (typically y = 0.005) a Born topology (0/0) 
event is generated. The corresponding cross section contains all virtual (V) corrections and is 
integrated over the phase space of all soft (S) photons with < y . Exponentiation of initial 
state radiation is implemented by modifying the 0(a) partial cross sections and interference 
terms in such a way that the derivative of the VS cross-section with respect to y is recovered in 
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the y — > limit. For example in the s-channel photon exchange contribution with initial state 
radiation: 

dQidy 2irs y 

exponentiation is carried out by the replacement u 2 + t 2 + u' 2 + t' 2 — > f (u 2 + 1 2 ) + u' 2 + 1' 2 
where / = 2 C v y^ e — 1. Events with hard photons are generated according to distributions 
where the soft photon eikonal factors are corrected by the Gribov-Lipatov [84] kernels. The 
relative probabilities of different topologies of final state photons are chosen according to the 
Poisson distribution: P(n\N) where n = n 7 — 1 and 

JV' = /3 e ln(l/y ), N F = ^HVvo) (7) 

A short description of program together with comparisons with other muon pair and wide angle 
Bhabha codes has been published [67]. A long write-up is also available [68]. 



OIMZ 


Z mass (GeV) 


OIMT 


Top quark mass (GeV) 


OIMH 


Higgs boson mass (GeV) 


OMAS 


a s (M z ) 


IOCH 


= 0(fi+fi-),= l(e+e-) 


IOEXP 


= 1 exponentiated , = 0(a) 


OW 


collision energy (GeV) 


0CTC1 


lower cos 6\+ in the ODLR frame 


0CTC2 


lower cos Q\- in the ODLR frame 


IOXI 


initial random number 



Table 22: Variables of the labelled common ICOM. 0CTC1,0CTC2 are used in setting up the LUT of the 
lepton scattering angles. To allow for the effects of the Lorentz boost the angular range should be chosen 
somewhat wider than that defined by the cuts in the LAB system. 

FEATURES OF THE PROGRAM: 

The execution of the program has three distinct phases: initialisation, generation and termina- 
tion. In the initialisation phase all relevant electroweak quantities are calculated from the input 
parameters Mz, M t Mjj and a s . Also a number of look up tables for quantities such as the 
lepton scattering angle and photon energies are created for use in the subsequent generation 
phase. This process is relatively time consuming, so the user should not be surprised if there 
is some delay between the execution of the program and the start of event generation. In the 
generation phase events with unit weight are generated by the weight throwing technique. The 
corresponding 4- vectors are stored in common C4VEC. The user may apply arbitary cuts and 
produce weighted histograms in subroutine FUSER. Histograms of unit weight events may be 
produced in subroutine FHIST. In the final, termination, phase the input parameters are printed 
out together with the exact cross section a CUT and its error, together with all histograms and 
plots. 

HOW TO USE THE PROGRAM: 

The program has a very short main program containing definitions of the most important 
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u 2 + t 2 + u' 2 + t' 2 



a'2 



(6) 
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NPAR(1 
NPAR(2 
NPAR(3 
NPAR(4 
NPAR(6 
XPAR(1 
XPAR(2 
XPAR(3 
XPAR(4 
XPAR(9 
XPAR(10 
YMA 
YMI 
WTMAX 



1_, weak loop corrections ON, OFF 
2,3 parameterisations of had. vac. pol. 

0. 1.2 two-loop OLOL s m 2 t correction 
1,0 weak box diagrams ON, OFF 

1, two-loop terms oc m\ ON,OFF 
initial lepton charge (-1.D0) 

final lepton charge (-1.D0) 

final lepton colour (1.D0) 

final lepton mass (GeV) 

QCD correction to (non-6 quarks) 

QCD correction to 

maximum value of E^/Eb eam (0.99D0) 
minimum value of E 1 / Eb eam (0.005D0) 
maximum value of the event weight (2.0D) 



Table 23: Control parameters defined in SUBROUTINE BHAGENE3. Default values are underlined or 
given in parentheses. 

input parameters, which are stored in the labelled common block ICOM. These variables are 
described in Table 22 The execution of the program has three distinct phases: (i) Initialisation, 
(ii) Generation of a single unit weight event, (iii) Termination. Each of these phases is entered 
via a call to subroutine BHAGENE3 in the main program: 
CALL BHAGENE3(MODE,CTP1,CTP2,CTM1,CTM2,CTAC,EPO,EMO) 

MODE is set to —1,0,1 for the initialisation, generation and termination phases respectively. 
The other parameters of BHAGENE3 define the kinematical cuts to be applied to the generated 
events: 

CTP1 = minimum value of cos 

CTP2 = maximum value of cos 

CTM1 = minimum value of cos Q\- 

CTM2 = maximum value of cos Q\- 

CTAC = maximum value of cos (j) co i 

EPO = minimum energy of l + (GeV) 

EMO = minimum energy of l~ (GeV) 
All these cuts are applied in the laboratory (incoming e + ,e~ centre of mass) system. The 
angle (j) co i is the collinearity angle between the l + and the l~ (CATC = -1 for a back-to-back 
configuration). In the calls of BHAGENE3 with MODE = 0,1 only this parameter need be 
specified. Other initialisation parameters of interest to users are defined in BHAGENE3 itself. 
A list of the most important of these can be found in Table 23. 

AVAILABILITY: 

From Compure Physics Communications Program Library, see http : //www . cs . qub . ac.uk/ cpc 
for more details. 
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4.3 BHLUMI 



Institute of Nuclear Physics, Krakow, ul. Kawiory 26a 
jadachOcernvm. cern. ch 

Institute of Computer Science, Jagellonian University, Krakow 
erichterOcernvm. cern . ch 

Department of Physics and Astron., University of Tennessee and SLAC 
bflwOslac . Stanford . edu 

CERN and Institute of Nuclear Physics, Krakow, ul. Kawiory 26a 
wasmOcernvm . cern . ch 

GENERAL DESCRIPTION: 

The program is a multiphoton Monte Carlo event generator for low angle Bhabha providing four- 
momenta of outgoing electron, positron and photons. The first 0(a 1 )yFS version was described 
in ref. [85]. The actual version 4. 02. a includes several types of the matrix elements. The most 
important (9(a 2 )y™s matrix elements (M.E.) is based on the Yennie-Frautschi-Suura (YFS) 
exponentiation. This M.E. includes exactly the photonic first order and second order leading- 
log corrections. In the (9(a 2 )y™s M.E. the other higher order and subleading contributions 
are included in the approximate form. The detailed description exists for the version 2.02 in 
ref. [30]. For the differences between the versions 2.02 and 4.02 the user has to consult ref. [6], 
the README file in the distribution directory and comments in the main program of the 
demonstration deck [44]. The only difference between versions 4.02 and 4.02a is correction to 
an important bug 95a. In order to correct it one has to replace j in eq. (3) in Ref. [6] with 

4 2) o] = (7p + 7«)ln A + - 7 - ^ - - 7 ln(l - ft) - - 7 ln(l - a 1 ). (8) 

We also provide patch to correct this in the sorce code of the versions 4.02, see AVAILABILITY 
below. This correction can affect the result of the program typicaly 0.05%, up to 0.08% for 
some event selections. 

FEATURES OF THE PROGRAM: 

BHLUMI consists in fact of the three separate event generators: BHLUM4, OLDBIS and 
LUMLOG, where OLDBIS is an improved version of the OLDBAB written by Berends and 
Kleiss at PETRA times, and LUMLOG is an event generator with the inclusive many photons 
emission strictly collinear to momenta of incomming/outcogoing fermions. M.E. of OLDBIS is 
limited to 0(a 1 ) and M.E. of LUMLOG includes exponentiated and non-exponentiated electron 
structure functions up to (9(o: 3 )ll. BHLUM4 includes four types of the exponentited M.E.: 
0(a 2 )[j S , C( al )AB S an d f° ur types of the non-exponentited M.E.: C(a 2 )A,B 5 C(« 1 )a,b where 
the cases A and B correspond to two kinds of ansatz employed for modeling the 0(a 2 L), second 
order NNL, contribution. The BHLUM4 program includes vacuum polarization, s-chanel 7 and 
Z exchange contributions, see ref. [59] in the approximation suitable for the low angle (below 
O.lrad.) scattering. The BHLUM4 does no include so called up-down interferences. However, 
OLDBIS does include them so it can be used to check how small they are. 



AUTHORS: 
S. Jadach 

E. Richter-Was 

B.F.L. Ward 

Z. Was 
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HOW DOES THE CODE WORK: 

The program ia a full scale Monte Carlo event generator. A single CALL BHLUMI produces one 
event, i.e. the list of the final state four- moment a of electron, positron and photons encoded 
in the common block. Depending on switch in the input parameters the program provides 
event with the variable weight WTMOD or with constant weight WTMOD=l. In the constant 
weight mode the calculation is done for M.E. of the 0(a 2 )^ FS type. In the variable weight mode 
WTMOD corresponds to the above M.E. but the user has also acces to all six types of the M.E. 
listed above (and even more) and may perform in a single run calculation for various types of 
the M.E. The choice of one of three sub-generators BHLUM4, OLDBIS or LUMLOG is decided 
through one of the input parameters. Program requires initialization before producing first 
MC event. There are many input parameters. The most important ones define the minimum 
and maximum angle (t chanel transfer). For weighted events it is possible to cover the angular 
range down to zero angle but the program is realy designed for "double tag" acceptance. It is 
possible to stop and restart the program from the next event in the series. The distribution 
directory incudes example demonstrating how to do it. 

DESCRIPTION OF THE OUTPUT: 

Program prints certain control output. The basic output of the program is the series of the 
Monte Carlo events and the user decides by himself which events are accepted or rejected ac- 
cording to his favourite selection criteria. The total cross section in nanobarns can be calculated 
for arbitrary cuts in a standard way 

" = "o-^ E W t (9) 

Accepted Events 

where the sum of the weights (variable or constant) is over all accepted events, N is total 
number of generated events and a is a reference (normalization) cross section in nanobarns 
provided by the program at the end of the MC generation. In the analogous standard way one 
may obtain any arbitrary distribution properly normalized. 

AVAILABILITY: 

The program is posted on WWW at http : / /hpjmiady . if j . edu.pl in the form of ".tar.gz" file 
together with all relevant papers and documentation in postscript. The version 4. 02. a which was 
used to produce all numerical results in this workshop consists of the version 4.02 described in 
ref. [6] and of the error patch posted in the same location http://hpjmiady.ifj.edu.pl. After 
workshop the equivalent version 4.03 will be released. The new version of BHLUMI will also 
contain new version of LUMLOG with the final state bremsstrahlung which was used in in the 
table/figure 19 and improved version of the BHLUMI matrix element without exponentiation 
which was used in this table/figure. 
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4.4 BHWIDE 



AUTHORS: 

Institute of Nuclear Physics, Krakow, ul. Kawiory 26a 
jadachOcernvm. cern. ch 

Dept. of Physics and Astron., Univ. of Tennessee 
placzek@hephp02 . phy s . utk . edu 

Dept. of Physics and Astron., Univ. of Tennessee and SLAC 
bflwOslac . Stanford . edu 

GENERAL DESCRIPTION: 

The program evaluates the large (wide) angle Bhabha cross section at LEP1/SLC and LEP2 
energies. The theoretical formulation is based on O(a) YFS exponentiation, with O(a) virtual 
(both weak and QED) corrections taken from ref. [60,86] as formulated in the program AL- 
IBABA. The YFS exponentiation is realized via Monte Carlo methods based on BHLUMI-type 
Monte Carlo algorithm, which is explained in refs. [30,85]. Thus, we achieved an event-by-event 
realization of our calculation in which arbitrary detector cuts are possible and in which infrared 
singularities are canceled to all orders in a. A detailed description of our work can be found in 
ref. [69]. 

FEATURES OF THE PROGRAM: 

The code is a full-fledged Monte Carlo event generator, so that the final particle four-momenta 
for the entire e + e~ + n~f final state are available for each event, which may be generated as 
a weighted or unweighted event, as the user finds more or less convenient accordingly. Thus, 
it is trivial to impose arbitrary detector cuts on the events. If the user wishes, he/she may 
also use the original BABAMC [82,83] type of pure weak corrections (there is a simple switch 
which accomplishes this). The expected accuracy of the program, when all tests are finished, 
is anticipated at ~ 0.2% in the Z-region and ~ 0.1% in the LEP2 regime. 

HOW DOES THE CODE WORK: 

The code works entirely analogous to the MC event generator BHLUMI 2.01 described in 
ref. [30] . A crude distribution consisting of the primitive Born level distribution and the most 
dominant part of the YFS form factors, which can be integrated analytically, is used to gener- 
ate a background population of events. The weight for these events is then computed by the 
standard rejection techniques involving the ratio of the complete distribution and the crude dis- 
tribution. As the user wishes, these weights may either be used directly with the events, which 
have the four-momenta of all final state particles available, or they may be accepted/rejected 
against a constant maximal weight WTMAX to produce unweighted events via again standard 
MC methods. Standard final statistics of the run are provided, such as statistical error analysis, 
total cross section, etc. 

DESCRIPTION OF THE OUTPUT: 

Program prints certain control output. The basic output of the program is the series of the 
Monte Carlo events. The total cross section in nanobarns can be calculated for arbitrary cuts 
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in the same standard way as for BHLUMI, i.e. user may imposed arbitrary experimental cuts 
by rejection. 

AVAILABILITY: 

The program can be obtained via e-mail from the authors. It will be posted soon on WWW at 
http : //enigma .phys .utk . edu as well as on anonymous ftp at enigma. phys .utk. edu in the 
form of ".tar.gz" file together with all relevant papers and documentation in postscript. It will 
also be available via anonymous ftp at enigma .phys .utk . edu in the directory /pub/BHWIDE. 



4.5 NLLBHA 



AUTHORS: 

Joint Institute for Nuclear Research, Dubna, 141980, Russia 
arbuzov@thsunl.jinr.dubna.su 

Joint Institute for Nuclear Research, Dubna, 141980, Russia 
kuraev@theor.jinrc.dubna.su 

CERN TH Division, Universita di Parma, INFN Sezione di Milano 
trenta@vxcern.cern.ch 

GENERAL DESCRIPTION: 

NLLBHA is a semi-analytical program for calculations of radiative QED and electroweak correc- 
tions to the small-angle Bhabha scattering at high energies. It takes into account complete 
(relevant at small angles) first order QED and electroweak corrections, the leading and next- 
to-leading QED corrections to 0(a 2 ) and the leading logarithmic contributions to 0(a 3 ). The 
corrections due to photon emission as well as the ones due to pair production are included. The 
theoretical uncertainty of the calculations is less then 0.1%. 

FEATURES OF THE PROGRAM: 

NLLBHA integrates numerically analytical formulae [2,8,23]. For the Born cross-section an ex- 
pansion for small scattering angles is used. The contributions due to real particle emission 
are integrated over symmetrical detector apertures. For the case of asymmetrical detectors 
(narrow-wide case) leading logarithmic contributions are calculated (next-to-leading are es- 
timated to be equal or less the ones in the narrow-narrow case). Cuts on the final particles 
energies are possible. Calorimetric set-up as well as other special experimental conditions are 
not implemented. 

HOW DOES THE CODE WORK: 

The code consists of the main part and of a series of subroutines which calculate separately 
radiative correction (RC) contributions from different Feynman diagrams and configurations. 
In the main part the flags, the parameters and the constants are defined. Using the flags one 
can define with their help the event selection (BARE1 symmetric or asymmetric are possible 
only), the order of corrections, switch on or off different contributions (like Z-boson exchange, 
vacuum polarization and light pair production). Then the user have to set the parameters: the 
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beam energy, the angular range, the energy cut. The electroweak parameters are calculated 
with the help of the DIZET package [87]. 

DESCRIPTION OF THE OUTPUT: 

At firdt the code prints the information about the chosen set-up, vacuum polarization (on/off), 
Z-boson contribution (on/off). Then the code prints the beam energy, the angular range and 
the electroweak parameters. After calculations it prints, for each value of x c (energy cut), the 
Born and the radiatively corrected (to different orders and approximations) cross-sections in 
[nb]; It also prints a line with the values of the different corrections in percent with respect to 
the Born cross-section. The normalizations and definitions used do directly correspond to the 
ones given in [2] where also the origin of all RC contributions can be found. 

AVAILABILITY: 

The code is available upon request from the authors. 

4.6 SABSPV 

AUTHORS: 

M. Cacciari DESY, Hamburg, Germany 

cacciariOdesy . de 
G. Montagna University of Pavia, Italy 

montagna@pavia.pv . infn . it 
O. Nicrosini CERN - TH Division (Permanent address: INFN Pavia, Italy) 

nicrosiniOvxcern . cern. ch, nicrosini@pavia.pv . infn . it 
F. Piccinini INFN Pavia, Italy 

piccininiOpavia . pv . infn . it 

GENERAL DESCRIPTION: 

SABSPV evaluates small angle Bhabha cross sections, in the angular region used for lumi- 
nosity measurement at LEP, and large angle Bhabha cross sections at LEP2. The theoretical 
formulation is based on a suitable matching between an exact fixed order calculation and the 
resummation of leading log radiative effects provided by the structure function techniques. 
The matching between the all-orders leading-log cross section, (T^) , given by the convolution 
of structure functions with kernel (Born) cross sections, and the O(a) one is realized according 
to the following general recipe: the order-a content of the leading-log cross section is extracted 
by employing the O(ct) expansions of the structure functions, thereby yielding cr£^. Denoting 
by cr s+v (k ) the cross section including virtual corrections plus soft photons of energy up to 
E 1 = k E, and by cr H (k ) the radiative O(a) cross section, the fully corrected cross section can 
finally be written as 

°a = 4l } " 4l } + * S+V (ko) + a H (k ) . (10) 

Equation (10) is in the additive form. A factorized form can also be supplied. It has the same 
O(a) content but also leads to the so-called classical limit, according to which the cross section 
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must vanish in the absence of photonic radiation. It reads 



a F = (1 + OnlJ^ll > ^nl = = > I 11 ) 



a being the Born cross section; C§ L contains the non-log part of the O(a) cross section, rep- 
resented by crjy^- 

In order to be flexible with respect to the different kinds of experimental cuts and triggering 
conditions, it makes use of a multi- dimensional Monte Carlo integration with importance sam- 
pling. A detailed description of the formalism adopted and the physical ideas behind it can be 
found refs. [2,46] and references therein. 

FEATURES OF THE PROGRAM: 

The code is a Monte Carlo integrator for weighted events. At every step, two kinds of events 
are fully accessible: 

(A) "two-body" events: they include tree-level events and radiative events in the collinear 
approximation; in this last case, information concerning the equivalent photons lost in the 
beam pipe is available, 

(B) three-body events: they include the radiative events e + e~ — > e + e~7 beyond the collinear 
approximation. 

No explicit photons beyond O(a) are generated; on the generated events, every kind of cuts can 
be imposed. O(a) corrections are available for the ^y(t)^(t) contribution (see for instance [72] 
for the soft plus virtual corrections and [88,89] for the hard bremsstrahlung contribution); 
all the other channels are treated in the leading logarithmic approximation 9 . This theoretical 
framework does exploit the fact that the ^y(t)^(t) channel is by far the most dominant one. It 
is therefore sufficient to evaluate exact order a corrections for this channel only. The other 
channels, which at the Born level contribute at the level of one per cent in the small angle 
region and of some per cents in the large angle region at LEP2 energies, can be evaluated in the 
leading log approximation. Higher order corrections are implemented in the structure function 
formalism [2]. The overall accuracy of the predictions performed by the code is, generically, of 
the order of 0.1% in the small angle regime and of the order of 1% in the large angle regime at 
LEP2 energies. 

HOW DOES THE CODE WORK: 

The code generates random integration variables within the "fiducial" cuts supplied via the 
input card (see below). These values are passed to the kinematics subroutines, which construct 
the full quadrimomenta for electron, positron and photon. The quadrimomenta are then fed 
to a trigger routine, which either accepts or rejects the event according to the cuts specified in 
it by the user. The control is then returned to the main integration routine, which generates 
weighted events, accumulates the cross section result for each single contribution and compose 
them as described in eqs. (10) and (11). Once in a given number of events the integrations 



9 Actually, in the present version of the program the up-down interference contribution is neglected. This is 
of no practical relevance for the small angle cross section, whereas it introduces an error of the order of some 
per mil in the large angle cross section at LEP2. 
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results and the related error estimates are evaluated and written to the output file. The error is 
also compared to the accuracy limit required, and the run stops when the latter is reached. The 
program can be restarted from its own output file, by specifying the same "physical" inputs 
and either a larger number of events or a higher accuracy. 

INPUT CARD: 

The following data card has to be provided via standard input: 



46.15D0 

24.D-3 58.D-3 0.D0 

24.D-3 58.D-3 0.D0 

0. 5D0 l.D-2 0.D0 0.D0 0.D0 
1 

1 

1. D5 0.D0 'SABSPV.OUT' 



EBEAM 

T1MIN, T1MAX, E1MIN 
T2MIN, T2MAX , E2MIM 
CAL0INPUT(1. . .5) 
ISIM 
ICALO 

EVTS, ACCLIM, IRESTART, OUTFILE 



These parameters have the following meaning: 

46.15D0 - the electron and positron beam energy, EBEAM. 
[~2~[ 24 . D-3 , 58 . D-3 , . DO - the electron minimum and maximum scattering angle (in radians) 
and the minimum electron energy (in GeV), T1MIN, T1MAX, E1MIM. These cuts are to be 
interpreted as "fiducial" cuts within which the events are generated, before going through the 
triggering routine. 



3j the same for the positron, T2MIM, T2MAX , E2MIM. 

4 0.5D0, l.D-2, 0.D0, 0.D0, 0.D0 - inputs that may be required by the cutting routines 
for the triggers. These values are stored in the vector CAL0INPUT(5) via the common block 
COMMON/CALOS. 

1 - flag for symmetric cuts, ISIM. The user has to specify if the experimental cuts asked 



for are (l) or not (0) symmetric for electron-positron exchange. If they are, choosing 1 saves 
computing time. 

|~6~| 1 - flag for choosing the triggering routine, ICALO. 

\Y\ 1.D5, 0.D0, 0, ' SABSPV.OUT' - these are inputs related to the Monte Carlo integration 
and to the management of the output. Namely, the total number of events, EVTS, the relative 
accuracy limit aimed at, ACCLIM, the restarting flag, IRESTART (if 1 the program tries to restart 
execution from the indicated output file, if it reinitializes it), and the output file name, 
OUTFILE. 

DESCRIPTION OF THE OUTPUT: 

The output file OUTFILE contains a description of the inputs provided to the code, the results 
of the Monte Carlo integrations for the various contributions and the final results with their 
standard statistical error. Moreover informations concerning the random number generator and 
the cumulants, that can be used to restart the program from where it stopped, are provided. 

AVAILABILITY: 

The code is available upon request to one of the authors. 
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4.7 UNIBAB 



AUTHORS: 

H. Anlauf TH Darmstadt & Universitat Siegen, Germany 

anlauf Scrunch. ikp .physik. th-darmstadt . de 
T. Ohl TH Darmstadt, Germany 

ohlOcrunch . ikp . physik . th-darmstadt . de 

GENERAL DESCRIPTION: 

UNIBAB is a Monte Carlo event generator designed for large angle Bhabha scattering at LEP and 
SLC energies. In its original incarnation [90,91], it was a simple QED dresser describing only 
multiphoton initial-state radiation, thus focusing on the exponentiation of soft photons and the 
resummation to all orders of the leading logarithmic corrections of the form (a/"7r) n \n n (s/ml). 
The first published version, UNIBAB version 2.0 [70] contains improvements in the exclusive 
photon shower algorithm used for the description of initial-state radiation, and many enhance- 
ments, such as final-state radiation using a similar photon shower algorithm. An electroweak 
library based on ALIBABA [60] was added. Initial and final state corrections are implemented 
in a fully factorized form. Version 2.1 of the program features the inclusion of longitudinal 
beam polarization. During this workshop the current version 2.2 was developed, which uses an 
implementation of the final state photon shower based on the exact lowest order matrix element 
for the process Z — > f f~f. Also, the electroweak library has been updated slightly to include 
the leading ra^-dependence and higher order QCD corrections to the Z width, as discussed in 
detail in [1]. 

FEATURES OF THE PROGRAM: 

The event generator UNIBAB calculates the QED radiative corrections through a photon shower 
algorithm. The actual implementation is based on an iterative numerical solution of an Altarelli- 
Parisi type evolution equation for the electron structure function. The effective matrix element 
for photon emission from the initial state assumes a factorized form of the radiative matrix 
element. Therefore it is exact for collinear emission. It also allows to generate finite transverse 
momenta of the radiated photons. For final state radiation, the algorithm employs an iterated 
form of the first order matrix element for Z — > f f~f, which gives a reasonable description 
of exclusive distributions that are sensitive to the details of the approximations used for the 
multiphoton matrix element, such as acollinearity distributions on the Z peak. 
UNIBAB generates only unweighted events. It is implicitly assumed that all scales in the hard 
subprocess are of the same order of magnitude, and the program does not yet include initial- 
final interference, thus the program is generally limited to the large angle region. Numerically, 
the effects from initial-final interference are sufficiently small in the vicinity of the Z peak. For 
details see the long write-up [70]. 

HOW DOES THE CODE WORK: 

UNIBAB consists of two layers, an external layer with a very simple user interface that allows easy 
interactive and batch control of the program, and an internal layer with a low level interface 
to the internal routines. It is however recommended to use the high level interface which 
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automatically takes care of parameter dependencies and properly reinitializes the Monte Carlo 
when a physics parameter is modified. 

In order to run the program, one has to specify several steering parameters that are internally 
translated into Monte Carlo parameters. The actual physical cuts have to be implemented in 
an external analyzer. The essential steering parameters are: 



• ctsmin, ctsmax: cuts on cos where 8* is the scattering angle in the boosted subsystem 
after taking initial state radiation into account. 

• ecut: minimum energy of the final state fermions. 

• acocut: maximum acollinearity of the outgoing e + e~ pair. 



An interactive run may look like: 

set ebeam 45.65 # Beam energy in GeV 

set masslz 91.1888 # Z mass 

set mass It 174 # top quark mass 

set masslh 300 # Higgs mass 

set ctsmin -0.8 

set ctsmax 0.8 

set ecut 20 

set acocut 30 # acollinearity cut in degrees 

init 

generate 100000 

close 

quit 

Additional switches control the inclusion or omission of certain contributions like weak box 
diagrams or ^-channel diagrams. For more details please consult the manual. 

DESCRIPTION OF THE OUTPUT: 

UMIBAB stores the generated events and all supplementary information for analysis (cross sec- 
tion, Monte Carlo error) in the proposed standard /hepevt/ common block [72] and must be 
read from there by a suitable analyzer. A simple yet very flexible tool for implementing a "the- 
orist's detector" is given by HEPAWK [92,93], which easily allows to obtain arbitrary distributions 
from the generated events. 

AVAILABILITY: 

The current version of UMIBAB may be downloaded via anonymous ftp from 
ftp : / / crunch. ikp . physik. th-darmstadt . de/pub/ anlauf /unibab 

along with up-to-date documentation. At the time of this writing (and for historical reasons), 
the program source and accompanying files are still distributed in the CERN patchy format. 
Platform-dependent Fortran77 source files will be made available upon request. For the sample 
test run, UMIBAB has also to be linked with the analyzer HEPAWK [92,93]. A more modern (auto- 
configuring) and self-contained version of the Monte Carlo generator will be made available in 
a future release after the end of the workshop. 
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5 Conclusions and outlook 



In this WG, the first systematic comparison of all the existing Monte Carlo event generators 
for the Bhabha process at LEP1 and LEP2 has been performed. This is one of our main 
achievements. The other one is that, as a result of these comparisons, the theoretical error 
of the small-angle Bhabha process is now reduced from 0.16% to 0.11% for typical LEP1 
experimental ES's, at the angular range of 1° — 3°. In parallel, an estimate of the theoretical 
error of the small-angle Bhabha process at LEP2 has also been fixed at 0.25%, for all possible 
experimental situations. The theoretical precision of the small-angle Bhabha scattering should 
be still improved by a factor of two at LEP1, in order to match the experimental precision. 
From the analysis performed, we conclude that a theoretical error of the order of 0.06% is 
reasonably feasible at LEP1, and the present study offers a solid ground for the next step in 
this direction. 

As far as the large-angle Bhabha process is concerned, the main result of this WG is that 
now we have comparisons not only among the semi- analytical benchmarks ALIBABA and 
TOPAZ0, but also among Monte Carlo event generators and on the Monte Carlo codes versus 
semianalytical programs. In spite of the fact that the comparisons involving Monte Carlo's 
do not change the conclusions of the previous LEP1 WG on the theoretical precision of large- 
angle Bhabha at LEP1 (see [1]), they give information about the performances of the Monte 
Carlo event generators themselves. In particular (except for some programs which have to be 
improved, either on the QED libraries or on the pure weak ones), the situation at LEP1 is 
generally under control with respect to the present experimental accuracy, both on and off Z 
peak. As far as LEP2 is concerned, a general agreement of the order of 2% has been achieved. 
There is certainly room for further improvements on this item, but for practical purposes the 
situation can be considered satisfactory. 
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